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Abstract 


Development  and  use  of  new  materials  and  technologies  for  use  in 
high  temperature  environments  are  dependent  on  the  availability  of  a 
material  model  appropriate  for  such  applications.  The  Bodner-Partom 
constitutive  model  shows  promising  flexibility  in  modeling  materials 
under  such  conditions.  Evaluation  of  the  model's  adequacy  and  range  of 
application  has  thus  far  been  hampered  by  a  lack  of  physical  understand¬ 
ing  of  the  model's  coefficients  and  the  lack  of  a  systematic  method  for 
evaluating  them. 

/ 

1/  n 

This  study  uses  experimental  data  for  IK1 718  at  1200  F  to  establish 
a  systematic  approach  to  the  evaluation  of  the  Bodner  coefficients. 
Parametric  studies  were  used  to  develop  a  comprehensive  description  of 
each  coefficient  and  its  impact  on  the  material  model.  A  simple  exten¬ 
sion  of  the  Bodner  constitutive  equations  to  model  cyclic  loading  situ¬ 
ations  was  carried  out. 

Results  show  that  systematic  evaluation  of  the  Bodner  coefficients 
from  experimental  data  is  achievable  although  somewhat  subjective  in 
nature.  Cyclic  load  studies  revealed  that  the  Bodner  model  is  capable 
of  modeling  isotropy  and  directional  strain  hardening  as  well  as  strain 
sof  tening. 


I.  Introduction 


The  continued  and  growing  use  of  high  performance  gas  turbine 
engines  on  military  aircraft  has  placed  the  United  States  Air  Force  in 
a  position  where  more  information  is  needed  on  material  behavior  under 
the  conditions  of  these  applications  [lj.  High  temperature  engine  oper 
ation  is  typified  by  complex  material  behavior  where  time  dependence 
and  temperature  effects  interact  with  a  varying  load  spectra  [2,3].  An 
analytical  technique  to  model  material  behavior  under  these  conditions 
is  required  to  allow  further  development  and  use  of  materials  and 
technology. 

Development  of  an  appropriate  analytical  model  has  been  a  slow 
process.  A  model  is  proposed;  then  it  must  be  evaluated,  compared  to 
experimental  data,  and  revised.  As  technology  advances,  new  materials 
and  new  applications  require  changes  in  the  models.  Today's  constitu¬ 
tive  models  have  evolved  from  classical  plasticity  over  a  period  of 
many  years. 

Early  material  models  were  based  on  the  assumptions  of  classical 
plasticity.  These  models  are  time  independent.  They  are  thus  unable 
to  adequately  describe  high  temperature  material  behavior  where  time 
and  rate  dependence  are  significant. 

The  classical  elasto-viscoplastic  models  followed.  They  corrected 
some,  but  not  all,  of  the  deficiencies  of  the  classical  plasticity 
models.  Classical  elasto-viscoplastic  theories  allow  for  strain  rate 
dependent  plastic  flow,  but  rely  on  a  strain  rate  independent  yield 
criterion  [2l]  (e.g.  von  Mises) .  Thus,  they  too  inadequately  describe 


high  temperature  material  behavior. 

Several  non-classical  elasto-viscoplastic  constitutive  models  have 
been  proposed  to  allow  analysis  of  the  behavior  of  rate  sensitive 
materials.  Each  attempts  to  model  the  mechanism  of  rate  dependence, 
temperature  dependence,  creep  and  relaxation  behavior,  and  each  predicts 
a  slightly  different  material  response  because  the  mechanisms  are 
modeled  differently.  A  brief  look  at  some  of  the  models  serves  to  point 
out  some  of  the  differences  between  them.  This  is  by  no  means  a  com¬ 
prehensive  overview.  It  is  merely  intended  to  provide  a  rudimentary 
understanding  of  some  of  the  model  types  in  use  today. 

Cernocky  and  Krempl  [24]  proposed  a  constitutive  model  that  relies 
on  an  overstress  mechanism  (a  variation  on  that  proposed  by  Malvern  in 
1951  [27])  and  an  equilibrium  material  stress-strain  response.  The 
difference  between  the  flow  stress  and  the  equilibrium  stress  is  called 
the  overstress.  The  model  is  non-linear  in  stress  and  strain  and  linear 
in  stress  and  strain  rates.  It  does  not  use  a  yield  surface.  This 
model  incorporates  stress-rate  history  effects  but  not  strain-rate 
history  effects  [17],  and  is  thus  unable  to  model  observed  material 
behavior. 

Lin  and  Wu  [23]  developed  constitutive  equations  based  on  the  endo- 
chronic  theory  attributed  by  Valanis.  In  the  endochronic  theories  of 
viscoplasticity,  a  transformation  is  used  to  move  from  the  absolute  time 
measured  by  a  clock  to  a  subjective  time  scale  that  is  essentially  a 
material  property.  Thus  the  current  value  of  stress  is  determined  by 
"the  memory  of  the  material  with  respect  to  the  endochronic  time"  [23]. 
Endochronic  theories  are  most  often  used  in  the  study  of  brittle 
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materials,  such  as  concrete.  This  model  does  not  rely  on  a  yield  sur¬ 
face.  It  models  creep  behavior,  but  does  not  predict  relaxation  be¬ 
havior  [17 ] ;  thus  it  is  unable  to  model  observed  material  response. 

The  model  proposed  by  Wu  and  Yip  [25]  is  also  an  endochronic  theory 
Unlike  that  of  Lin  and  Wu,  this  model  predicts  both  creep  and  relaxation 
behavior.  However,  where  most  endochronic  theories  do  not  require  use 
of  a  yield  surface,  this  one  does;  the  internal  measure  of  time  depends 
on  the  plastic  strain  rate.  Thus,  use  of  this  model  requires  that  a 
yield  criterion  be  established. 

The  Bodner-Partom  model  for  plastic  flow  [22]  does  not  depend  on  a 
yield  surface.  It  predicts  creep  and  relaxation  behavior;  it  incorpor¬ 
ates  time,  rate  and  temperature  effects.  Both  stress  and  strain  rate 
history  effects  are  modeled.  And  it  can  be  used  to  model  material 
behavior  during  application  of  cyclic  loads  as  well  as  monotonic  loading 
conditions. 

Because  the  Bodner-Partom  flow  law  appears  to  incorporate  mechan¬ 
isms  for  modeling  these  important  parameters  of  material  behavior,  it 
has  a  great  potential  as  an  analytical  tool.  At  this  time,  however, 
there  is  no  established  method  for  determining  the  eight  material  co¬ 
efficients  used  in  fffcfis  model.  Efforts  to  evaluate  the  adequacy  and 
flexibility  of  the  model  and  to  use  it  effectively  have  been  hampered 
by  the  non-availability  of  the  required  material  parameters. 

D.  C.  Stouffer  [9]  evaluated  the  coefficients  used  in  the  Bodner 
flow,  law  to  model  Inconel  100.  The  data  base  he  used  for  this  work  con¬ 
sisted  of  twenty  experiments,  including  both  tensile  and  creep  tests. 
Although  he  did  develop  appropriate  values  for  the  coefficients,  he  did 


so  without  developing  much  physical  understanding  of  them.  It  was  not 
known  how  sensitive  the  model  may  be  to  small  variations  in  the  values 
of  the  various  coefficients.  There  was  little  guidance  as  to  what  data 
should  be  used  to  evaluate  which  coefficients,  and  no  guidance  as  to 
how  sensitive  the  model  is  to  small  deviations  from  the  "best  fit" 
values  at  any  given  point. 

This  thesis  provides  a  systematic  method  for  evaluating  the  coeffi 
cients  of  the  Bodner  model  from  experimental  data.  The  work  done  by 
Stouffer  [9]  is  used  as  a  starting  point;  some  guidance  is  given  regard 
ing  what  constitutes  an  adequate  data  base  for  the  evaluation  process. 

A  comprehensive  discussion  of  the  coefficients  and  their  characteristic 
effects  on  the  Bodner  model  provides  additional  insight  into  the  model. 
A  study  of  the  Bodner  model  under  conditions  of  cyclic  loading  esta¬ 
blishes  the  validity  of  such  uses  of  this  model  and  suggests  additional 
characteristics  that  must  be  incorporated  into  the  cyclic  formulation 
to  allow  adequate  representation  of  material  behavior  under  such  condi¬ 


tions. 


in  this  formulation  that  the  plastic  deformations  are  described  by  the 
functional  relationship 


where  is  written  explicitly  as  a  function  of  the  deviatoric  stress 

J~,  absolute  temperature  T  and  a  set  of  internal  state  variables  Z,  . 

k 

The  current  values  of  the  Z^  are  meant  to  reflect  all  the  pertinent 
effects  of  any  prior  inelastic  deformations.  Thus,  we  can  write 

X  -  (dJ  (J2,  T,  Zk)/J2y*  (5) 

The  form  chosen  by  Bodner  and  Partom  [22]  for  expressing 
(J2,  T,  Z^) ,  motivated  by  the  equations  of  dislocation  dynamics,  is 

D2  =  D02e*P  -[4]D  ^ 

where  Dq  is  a  limiting  value  for  the  plastic  strain  rate  in  shear  and  n 
is  a  temperature  dependent  state  variable  which  describes  strain  rate 
sensitivity.  The  variable  Z  is  a  material  hardness  parameter  indicating 
in  some  sense  the  material's  internal  resistance  to  plastic  flow,  and  is 
discussed  in  greater  detail  subsequently. 

Some  insight  into  the  form  of  Equation  (6)  may  be  gained  from  a 
brief  look  at  dislocation  dynamics.  Dislocation  velocity  cannot  in¬ 
crease  without  limit  as  stress  increases;  it  must  approach  some  limiting 
value.  Quite  generally  the  velocity  may  be  written 

V  -  v*P  (a  )  (7) 

m  s 

* 

where  v  is  the  dislocation  velocity,  v  the  limiting  value  of  dislocation 

velocity,  CTg  the  shear  stress  between  slip  planes  and  P^  the  average 

* 

probability  of  v  taking  on  the  value  v  at  a  given  instant.  Whatever 


functional  form  is  chosen  for  P  (a  ) ,  it  must  satisfy  two  conditions: 

ms 

P  (a  )  must  be  zero  .at  a  =  0,  and  P  (o  )  must  approach  unity  asymptoti- 
ms  .s  ms 

cally  for  high  values  of  ag.  The  form  proposed  by  Gilman  [13]  which 
satisfies  both  these  conditions  is 


_  -D/a 
P  =  e  s 
m 


where  D  is  a  drag  stress  constant,  sometimes  called  the  dislocation  drag 
term.  Note  the  similarity  between  the  resulting  expression 


*  r  D  ' 

=  v  exp  -  — 

O’ 

s  _ 


and  Equation  (6). 


If  Equations  (2),  (3),  (5),  and  (6)  are  combined,  the  resulting 
expressions  are  found,  considering  a  uniaxial  stress  a  : 


*  AS  *=  Xo 

XXX 


a  2 
x  D  exp 

^  0 
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The  hardness  parameter  Z  consists  of  two  terms,  as  does  the  time 
rate  of  hardening,  Z.  The  first  term  represents  hardening  due' to  plastic 
work,  W  ,  while  the  second  incorporates  the  effects  of  thermal  recovery 
of  the  hardening  generated  by  plastic  work.  By  examining  the  mechanics 
of  strain  hardening  and  recovery  on  a  microscopic  scale,  some  insight 
into  the  forms  of  these  two  terms  can  be  found. 

The  softer  a  material  is,  the  easier  it  is  for  dislocations  to 


propagate  through  the  material  lattice.  It  is  the  result  of  this  propa¬ 
gation  of  dislocations  that  is  called  "plastic  strain."  Dislocations 
scattered  through  the  material  are  one  of  many  mechanisms  that  serve  to 
impede  the  propagation  of  mobile  dislocations.  Strain  hardening  is  the 
multiplication  of  dislocations  in  the  material  in  such  a  vay  as  to  Impede 
dislocation  motion. 

A  material  containing  these  "impedance"  dislocations  has  a  higher 
internal  energy  than  the  same  material  without  the  dislocations.  The 
higher  the  dislocation  density,  the  higher  the  internal  energy  [ll]. 

Over  time,  the  material  will  tend  toward  a  lower  internal  energy  state. 
The  material  reduces  its  dislocation  density  and  hence  its  internal 
energy  by  eliminating  or  relocating  impedance  dislocations.  This  reduc¬ 
tion  in  dislocation  density  is  what  is  called  recovery,  and  the  result 
is  a  softer  material. 

In  Bodner's  formulation,  [22],  the  work  hardening  expression  is 
written  in  terms  of  Z: 


(11) 


The  change  in  hardnes^  Z  with  respect  to  plastic  work  (per  unit 
volume)  is 


wh  «  m  (Z.  -  Z)  (12) 

3W  i 

P 


where  m  is  a  material  constant  describing  the  rate  of  work  hardening 
and  Z,  is  the  saturation  (maximum)  value  of  Z.  Note  also  that 


The  thermal  recovery  term  is  postulated 


Z 

rec 


=  -AZ, 


by  Bodner  [l]  to  have  the  form 
Z  -  Z„ 


(14) 


where  is  the  non-work  hardened  value  cf  Z  at  a  given  temperature,  and 
A  and  r  are  temperature  dependent  material  coefficients.  The  complete 
expression  for  the  rate  of  hardening  then  becomes 


Z  -  m  (Z  -  Z)  cep 


(15) 


Equations  (10)  and  (15)  together  are  assumed  to  provide  sufficient 
information  to  model  the  material  behavior.  For  test  conditions  where 
there  is  no  load  reversal,  a  forward  integration  technique  may  be  used 
to  incrementally  describe  material  behavior.  When  there  is  a  load  re¬ 
versal,  however,  proper  application  of  the  constitutive  equations  is 
considerably  more  complex. 

Under  cyclic  loading  the  hardness  parameter  Z  develops  directional 
characteristics.  Strain  hardening  caused  by  tensile  loading  affects  the 
material's  subsequent  behavior  under  compressive  loads.  Separate  values 

of  the  hardness  Z  would  thus  be  indicated  for  use  in  Equation  (10), 

•  *  • 

depending  on  the  sign  of  the  applied  stress.  When  the  applied  stress  is 
tensile,  the  applicable  Z  is  called  ZC;  when  the  applied  stress  is  com¬ 
pressive,  the  applicable  Z  is  called  Zc.  A  simple  model  [7]  incorpor- 
ating  isotropic  and  directional  hardening  allows  Z  and  Z  to  be  incre¬ 
mentally  computed  from  the  initial  value  Zq  according  to  a  few  basic 
principles. 


Initially, 


M 


Changes  from  this  state  are  computed  according  to 


Z  =  q  Z  +  (1  -  q)  Z  U 


Z  ■  q  Z  -  (1  -  q)  Z  U 


(17a) 


(17b) 


where  depends  on  the  sign  of  the  applied  stress,  and 

•  *p  ^eff  r 

z  =  m  (Z1  -  zeff)aep  -AZ1  — - 


where 


„t  _c 
z  =  Z  or  Z 
eff 


•  t  •  c  t 

depending  on  U^.  Using  these  computed  values  of  Z  and  Z  ,  values  of  Z 
and  Zc  are  determined  incrementally  over  time  using  the  relationships 


t  t  *  t 
Z  -  Z  +  Z  At 


Zc  *  Zc  +  ZCAt 


(20a) 


(20b) 


Thus,  at  the  end  of  the  first  incremental  time  step,  there  are  two 

a  t  • 

£  c 

values  of  hardness,  Z  and  Z  .  The  parameter  q  is  a  significant  feature 

of  this  formulation,  and  its  value  corresponds  to  the  uature  of  the 

hardening  process.  If  the  hardening  is  purely  isotropic,  q  ■  1;  for 
t  c 

this  case  Z  «Z  .  A  value  of  q  *  0  indicates  that  the  hardening  is  purely 
directional.  This  model  also  allows  for  situations  where  hardening  is 
pi.  ly  isotropic  and  partly  kinematic;  in  this  case  some  positive  frac¬ 
tional  value  is  assigned  to  q.  In  general  work  related  to  strain 


hardening  it  has  been  proposed  [7]  that  these  constitutive  equations 
can  also  be  associated  with  strain  softening;  in  this  case  a  value  of 
q<0  is  used.  In  the  subsequent  discussion  it  is  this  phenomenon  of 
strain  softening  that  will  be  dealt  with.  The  Bodner  constitutive  model 
is  capable  of  modeling  both  hardening  and  softening;  this  is  very  rare. 


Applications  for  Evaluating  Coefficients 


The  preceeding  paragraphs  briefly  present  a  technique  whereby 
Equations  (10)  and  (15)  can  be  used  to  describe  material  behavior  under 
cyclic  loading.  Under  certain  other  loading  conditions,  these  two 
equations  can  be  simplified  to  facilitate  the  evaluation  of  some  of  the 
coefficients. 

When  both  strain  rate  and  stress  are  constant,  as  is  the  case  when 
saturation  is  reached.  Equation  (10)  indicates  that  the  hardness  para¬ 
meter.  Z  must  also  be  constant.  Furthermore,  for  tests  of  short  duration 
(with  no  recovery)  the  material  must  be  in  a  fully  work  hardened  state 
at  saturation  "to  obtain  the  maximum  value  of  stress"  [9],  so  that 
Z  *  Z^.  For  these  steady  flow  conditions, 

J’-i  *8(?T] 

A?  r  n+l/ZlVnl 

~ar  ' exp  r  J 

-ln-2T  +  2 

ln[~ln_2rJ  J  J 

Inf- In  1  -  -2n  In  0  +  [2n  In  Z.  +  ln(~2n  jl 


(21) 


From  this  form  it  is  seen  that  In  (-lnC-^ —  ))  is  linearly  related  to  In 

o 


C.  Characterization  of  appropriately  plotted  experimental  data  will 
allow  determination  of  n  from  the  slope-  of  the  line.  Once  n  is  deter¬ 


mined,  can  be  deduced  mathematically. 


For  test  conditions  such  as  high  strain  rate  situations  which  pre¬ 
clude  recovery.  Equation  (15)  can  be  written 


dZ  =  m  (Z,  -  Z)  dW 
1  P 


(22) 


This  can  be  integrated  to  give 


In  (Z.  -  Z)  =  In  (Z,  -  Z.)  -  mW 

1  10  p 


(23) 


where  Z.  is  defined  as  the  initial  value  of  hardness  Z.  This  shows  that 


a  linear  relationship  exists  between  In  (Z^  -  Z)  and  the  plastic  work 


V  .  The  value  of  Z  for  each  data  point  of  interest  can  be  computed 
P 


from  an  inversion  of  Equation  (10)  which  reduces  to 


°E 


2D 

2n  .  o 
In 


n+1 


/3ep 


*5n 


(24) 


Plotting  In  (Z^  -  Z)  against  plastic  work  allows  the  determination  of 

a  value  for  Z^  from  jttie  linearly 'extrapolated  value  of  In  (Z^  -  Z)  at 

W  *  0.  The  slope  of  this  line  defines  a  value  for  m.  Once  the  initial 
P 

value  of  hardness  Z^  and  the  rate  of  work  hardening  m  are  known,  stress- 
strain  curves  can  be  generated  using  a  forward  integration  technique  on 
Equations  (10)  and  (15) .  This  allows  comparison  of  the  high  strain  rate 
analytical  prodictions  with  the  data. 

Based  on  the  values  determined  for  n,  m,  Z^  and  Z^,  stress-strain 
curves  are  generated  using  forward  integration.  This  procedure  is 


insensitive  to  the  values  chosen  for  A,  r,  and  (as  they  appear  only 
in  the  recovery  term)  so  the  fact  that  these  coefficients  have  not  yet 
been  determined  exactly  does  not  impact  these  curves.  Comparison  of 
these  curves  with  appropriate  experimental  data  allows  refinement  of  the 
values  chosen  for  n,  m,  and  Z^  to  best  model  experimental  data. 

A  parametric  study  of  the  coefficients  and  their  impact  on  the 
resulting  stress-strain  curves  provides  some  guidance  for  the  refinement 
process.  The  results  of  this  study  are  contained  in  Figures  1  through  4 
The  nominal  values  used  for  the  coefficients  are  given  in  Table  I: 
deviations  from  these  values  are  annotated  on  the  individual  figures. 


TABLE  I 

Seminal  Coefficient  Values  Used  in  Bodner  Equations 
For  Stress-Strain  Curves 


Variable 

Description 

Value 

n 

Strain  Rate  Sensitivity  Parameter 

3.0 

m 

Rate  of  Work  Hardening 

.219  MPa'1 

Z 

0 

Initial  Value  of  Material  Hardness  Z 

1669  MPa 

*i 

Saturation  (Maximum)  Value  of  Z 

1795  MPa 

Z2 

Minimum  Recoverable  Value  of  Z 

.  S 

718  MPa 

A 

Hardening  Recovery  Coefficient 

. 001  sec  1 

r 

Hardening  Recovery  Exponent 

7.0 

D 

o 

Limiting  Value,  Plastic  Strain  Rate 
in  Shear 

in6  "I 

10  sec 

-  E 

Elastic  Modulus 

162.5  x  10^  MPa 
1.3  x  10"5  sec-1 

£ 

Strain  Rate 

Figure  1  shows  the  impact  of  varying  n  and  at  a  strain  rate  of 
1.3  x  10  ^  sec  (This  strain  rate  is  one  of  those  for  which  experi¬ 

mental  data  exists).  The  four  stress-strain  curves  overlay  one  another 
quite  closely.  Any  n  value  in  the  range  from  r.  =  2.5  to  n  *  3.211 
clearly  predicts  material  behavior  quite  similar  to  that  predicted  by 
any  other  n  value  in  the  range. 

Figure  2  shows  the  impact  of  varying  m  plus  and  minus  ten  percent 
of  its  mean  experimental  value  (.219  MPa  ^) ,  again  at  a  strain  rate  of 
1.3  x  10  sec  Again,  the  four  curves  overlay  one  another  closely, 
but  this  time  the  effect  is  to  change  the'  shape  of  the  stress— strain 
curve  slightly  (increasing  m  causes  higher  saturation  stress  over  the 
domain  shown) . 

Figure  3  shows  the  impact  of  varying  Zq  plus  and  minus  ten  percent: 
of  its  mean  experimental  value  (1669  MPa),  at  a  strain  rate  of  1.3  x 
10  5  sec  1.  Changes  in  Zq  are  much  more  readily  seen  than  changes  in 
m  in  these  figures.  Increasing  ZQ  is  equivalent  to  increasing  the 
initial  hardness  of  the  material  without  changing  the  saturation  hard¬ 
ness  Z^.  Thus  the  nature  of  the  hardening  portion  of  the  curve  is 
affected.  Higher  values  of  Zq  correspond  to  a  flatter  appearance  of  the 
stress-strain  curve.-  5 

Figure  4  shows  the  impact  of  varying  m  and  Zq  at  a  strain  rate  of 
1.3  x  10  5  sec  The  shape  of  the  stress-strain  curve  can  clearly  be 
changed  significantly  by  different  choices  of  values  for  m  and  Zq. 

A  procedure  for  determining  appropriate  values  for  the  remaining 
material  coefficients  (A,  Z^  and  r)  requires  use  of  the  full  equations 
(equations  (10)  and  (15)).  To  demonstrate  more  vividly  the  character 


>  %  vvs 


Variation  of 


Variation  of  m  and  Z 


of  the  various  coefficients,  a  parametric  study  was  undertaken.  The 
information  contained  in  Equations  (10)  and  (15)  can  be  graphed  as  a 
non-dimensionalized  strain  rate  versus  a  non-dimensionalized  stress  to 
allow  for  a  more  general  presentation.  In  this  study,  nominal  values 
were  selected  for  all  coefficients  but  one,  and  various  values  for  this 
coefficient  of  interest  were  selected  to  illustrate  the  impact  of 
varying  that  parameter.  The  nominal' values  used  in  this  study  are 
presented  in  Table  II. 

TABLE  II 

Nominal  Coefficient  Values  Used  in  the  Bodner  Equations 
For  Strain-Rate  Versus  Stress  Curves 


Variable 

Description 

Value 

n 

Strain  Rate  Sensitivity  Parameter 

5.0 

m  times 

Dimensionless  Form  of  m 

1000. 

Vzi 

Dimensionless  Form  of  Z ^ 

.1 

A/D 

o 

Dimensionless  Form  of  A 

.000001 

r 

Hardening  Recovery  Exponent 

3.0 

Equation  (10)  tfas  non-dimensionalized  as  follows 


jiA  y  -2±i|  ,2a 
l 2  W  n  J  A? 

IK-)"  1  •  “B1 


where  a  is  defined  as  — .  For  this  study,  the  condition  of  saturation  was 
assumed,  meaning  that 

Z  *  0 


Substituting  this  into  Equation  (15)  leads  to  the  following: 


m  (Z1  -  Z)aep  -  A  Z 


Z  -  Z, 


Note  again  that  —  =  a.  Then  Equation  (26)  can  be  written 


m(Z1  -  aa)oep  =  AZ 


■  aa  -  Z2  1r 


mZ, 


«  *  T  oa  Z_  ]r 


Let  x  = 


-r  and  y 
L1 


ep .  Then 


Z„  .  r 


mZ^l-coOxy  *  ^  (ax -7^  ^ 

*  s  ,  0  1 

A  /  ^2  \r 

mZ1(l-ax)xy  -  —  (ax-cp  J 


(26) 


(27) 


Using  an  iterative  procedure,  a  value  for  x  (or  —  )  can  be  computed  for 

CP  1 

each,  value  of  y  (or  — ) .  Using  Equations  (25)  and  (27),  an  initial 

o 


estimate  of  the  proper  value  for  x  is  made.  Equation  (25)  is  substituted 


into  Equation  (27)  and  the  left  hand  side  of  Equation  (27)  is  calculated. 


As  this  will  not  be  zero  unless  the  value  for  x  is  guessed  exactly,  the 
difference  between  this  calculated  value  and  zero  serves  as  a  guide  to 
the  next  approximation  of  x.  At  high  strain  rates,  where  there  are  no 


significant  recovery  effects,  Z  =  Z^.  Under  these  conditions,  with  Z^ 
substituted  into  Equation  (25),  a  value  for  x  can  be  found  directly  from 


this  equation  alone. 

ep 

Graphing  the  dimensionless  strain  rate  —  versus  the  dimensionless 

Do 

stress  —  illustrates  the  impact  of  A,  r,  and  Z  on  the  predicted 

2 

material  behavior.  Figure  5  serves  as  a  guide  to  the  information  con¬ 
tained  in  this  type  of  figure;  it  is  based  on  the  studies  contained  in 
Appendix  A.  Figures  6,  7,  and  8  and  the  following  discussion  serve  to 
clarify  the  behavior  of  A,  r,  and  Z^. 

A/Dg  is  the  dimensionless  form  of  the  coefficient  in  th%  expression 
for  thermal  recovery  of  work  hardening  (Equation  (27)).  The  "branching" 
effect  in  the  lower  left  of  Figure  6  shows  that  the  primary  effect  of 
changes  in  A/D^  is  to  extend  the  predicted  recovery  regime  of  the  mate¬ 
rial.  The  onset  of  the  dominance  of  thermal  recovery  is  described  by 
this  variable.  For  larger  values  of  A/Dq,  thermal  recovery  dominates  at 
higher  strain  rates;  for  lower  values  of  A/Dg,  thermal  recovery  remains 
a  negligible  phenomefnbn  until  relatively  slower  strain  rates  occur. 
Changes  in  A/D^  have  no  significant  influence  at  strain  rates  faster 
than  that  at  which  recovery  becomes  significant.  As  a  part  of  the 
equation  descriptive  of  thermal  recovery,  A/Dq  does  not  effect  material 
behavior  at  higher  strain  rates  where  recovery  is  negligible. 

The  coefficient  r  appears  as  an  exponent  in  the  expression  for 


thermal  recovery  of  work  hardening  (Equation  (27)).  If  A  (or  A/D^) 


g  r 


HIGH  STRAIN  RATE  REGION  J/ 
WHERE  RECOVERY  EFFECTS  / 
MAY  BE  NEGLECTED  '  / 


SLOPE  DETERMINED 
BY  N 


MID  -  RANGE  : 
RECOVERY  AND 
HARDENING 


LOCATION  OF  INTER¬ 
SECTION  DETERMINED 
BY  A 


slope  Determined  by  r 


STRESS 


Figure  5.  Strain  Rate  Versus  Stress  Curve: 
A  Guide  to  the  Information 


were  to  be  written  a'fc’a  function  *of  the  onset  of  recoverability,  r  would 
be  the  curvature  of  that  function.  Where  the  value  of  A/Dq  determines 
at  which  strain  rate  recovery  begins  to  be  important,  r  determines  how 
quickly  recovery  becomes  a  dominant  characteristic  (see  Figure  7).  As 
with  A,  the  effect  of  changing  r  is  not  seen  at  strain  rates  too  high 
to  allow  thermal  recovery  of  plastic  work  to  be  significant. 


Z,/Z,  is  the  dimensionless  form  of  the  coefficient  Z_.  A  value  for 


.07  .28  49  .70  SI  1.12  133 

07z.| 

Figure  7,  Variation  of  r 

is  determined  mathematically  when  a  value  is  selected  for  n;  thus 
the  behavior  illustrated  in  Figure  7  is  that  of  Z^.  The  value  of  Z 2 
controls  where  the  curve  becomes  strain  rate  insensitive  at  very  low 
strain  rates.  The  value  of  Z ^  represents  the  minimum  recoverable  value 
of  hardness. 

A  more  complete  discussion  of  the  material  coefficients  is  to  be 


found  in  Appendix  A. 


Applications  to  IN  718 


Use  of  the  constitutive  equations  developed  in  the  preceeding 
section  to  model  material  behavior  requires  that  appropriate  values  be 
determined  for  the  coefficients.  In  this  exercise,  the  material  to  be 
modeled  was  Inconel  718,  at  1200°F.  The  data  was  provided  by  the 
Materials  Laboratory,  Air  Force  Wright  Aeronautical  Laboratories.  Con¬ 
stant  strain  rate  tensile  test,  creep  test  and  cyclic  loading  test 
data  were  provided. 

/3ep 

The  linear  relationship  between  In  (-In and  In  a  for  high 

o 

strain  rates  (Equation  (21))  provides  a  starting  point  for  the  deter¬ 
mination  of  the  material  coefficients  of  the  Bodner  model.  Equation 
(21)  is  repeated  here  for  continuity: 


Jr \  IP  nj.1 

ln(-ln(  2d"))  *  “2n  ln  O  +  (2n  In  +  In  — — ) 
o 


Considering  the  case  of  high  strain  rates,  where  recovery  effects  can  be 

ignored,  allows  n  and  Z^  to  be  evaluated.  If  the  appropriate  experi- 

/3ep  _ 

mental  data  is  graphed  as  In  (-ln(— ^— ))  versus  In  a,  and  a  straight 

o 

line  approximation  made  to  the  data,  the  slope  of  the  line  will  be  -2n. 

«  «  * 

A  value  for  can  then  be  determined  mathematically,  using  Equation  (21) 
and  the  coordinates  of  a  point  on  the  line.  Plotting  the  high,  strain 
rate  experimental  data  will  thus  allow  n  and  Z^  to  be  evaluated. 

The  full  data  base,  tabulated  in  Table  III,  was  graphed  in  Figure  9. 
This  allowed  an  initial  evaluation  to  be  made  regarding  which  data 
points  are  representative  of  situations  with  high  strain  rates  and  no 
recovery.  Specifically,  the  three  data  points  of  Figure  9  which  came 
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TABLE  III 


Experimental  Data  Base;  IN  718,  1200°F 


Stress  (KSI) 

Strain  rate 

Comments 

80 

_q 

4.0x10  sec 

Creep 

test 

100 

1.6x10  ^sec  ^ 

Creep 

test 

105 

1.8x10  ^sec‘^ 

Creep 

test 

110 

3.0x10  ^sec  ^ 

Creep 

test 

no 

3.0x10  ^sec  ^ 

Creep 

test 

no 

—8  —1 

3.0x10  sec 

Creep 

test 

no 

2.8x10  ^sec  ^ 

Creep 

test 

no 

3.6x10  ®sec  ^ 

Creep 

test 

120 

1.7x10  ^sec  ^ 

Creep 

test 

• 

120 

1.4x10  ^sec  * 

Creep 

test 

120 

1.5xl0~^sec  ^ 

Creep 

test 

125 

4.0x10  ^sec  ^ 

Creep 

test 

125 

3.0x10  ^sec  ^ 

Creep 

test 

130 

8.3x10  ^sec  ^ 

Creep 

test 

134 

1.1x10  ^sec  ^ 

Constant  strain 

rate 

135 

.  i 

5.5x10  ^sec  ^ 

Creep 

test 

140 

7.0x10  ^sec  ^ 

Creep 

test 

142 

1.3x10  ^sec  ^ 

Constant  strain 

rate 

147 

1.6x10  ^sec  ^ 

Constant  strain 

rate 

27 


The  two  remaining  data  points  were  graphed  on  In  (-ln(  ))  versus 

•  p 

In  a  axes,  where  e  was  considered  to  be  equal  to  the  total  constant 
strain  rate  £  in  each  case.  A  straight-line  was  drawn  to  characterize 
this  data.  The  results  are  depicted  in  Figures  10  and  57. 

According  to  Equation  (21) ,  the  slope  of  this  line  was  twice  the 
negative  value  of  the  coefficient  n.  Three  values  for  n  were  determined 
to  be  fairly  characteristic  of  the  data.  The  slope  was  first  determined 
using  the  coordinates  of  the  two  data  points,  resulting  in  a  value  of 


lN  CT  (MPA.) 


Figure  10.  Determination  of  n  and  Z 


n  =  3.211.  Working  graphically,  to  evaluate  the  effect  on  accuracy  of 
a  graphical  technique,  a  second  slope  was  measured  which  yielded  a  value 
of  n  =  3.125.  For  comparison,  the  line  representing  a  value  of  n  =  3.0 
was  drawn  on  the  graph  (see  Figure  57);  within  the  limits  on  accuracy 
imposed  by  the  graphical  methods,  this  line  also  seems  fairly  represen¬ 
tative  of  the  data.  For  a  final  point  of  comparison,  the  line  repre¬ 
senting  n  =  2.5  was  also  drawn.  This  line  is  not  so  different  from  the 
lines  representing  the  other  values.  It  would  appear  that  the  proper 
discrimination  of  the  best  value  for  n  is  not  possible  from  this  repre¬ 
sentation  of  the  data.  The  four  values  for  n  were  carried  forward  — 
two  limiting  values  (2.5  and  3.211)  and  the  two  intermediate  values 
(3.0  and  3.125)  —  in  the  hope  that  further  comparison  with  the  data 
would  provide  a  more  solid  basis  for  choosing  a  value  for  n.  - 

For  each  value  of  n,  a  corresponding  value  for  was  computed 
using  Equation  (21).  A  more  detailed  description  of  the  determination 
of  values  for  n  and  is  to  be  found  in  Appendix  B. 

The  second  study  undertaken,  to  determine  values  for  m  and  Z^,  was 
prompted  by  the  linear  relationship  between  In  (Z^  Z)  and  the  inelastic 
work  W  at  high  strain  rates  (Equation  (23),  repeated  here  for  conven- 

;  f 

ience) : 

ln(Zx-Z)  =  IoU^Zq)  -  mWp  .  (23) 

Again,  data  from  high  strain  rate  tensile  tests  was  used,  this  time  to 

evaluate  m  and  Z^.  The  data  used  in  this  step  came  from  a  single  con- 

•  -5  -1 

stant  strain  rate  test  at  £  =  1.3  x  10  sec  .  This  data  was  chosen 
because  the  strain  rate  was  fast  enough  to  allow  use  of  Equation  (23)  — 


which  neglects  recovery  effects  —  and  because  the  experimental  stress- 
strain  curve  was  well  defined  and  well-behaved  (see  Figure  11). 


Evaluation  of  the  data  began  with  the  computation  of  the  plastic 
strain  rate  indicated  by  the  data.  The  line  representing  Young's 
Modulus  (23.6  x  10  KSI)  [10]  was  drawn  on  the  experimental  stress- 
strain  curve  and  the  origin  translated  to  meet  the  initial  condition  of 
zero  strain  at  zero  stress.  The  values  for  total  strain  at  selected 
stress  levels  were  read  off  the  graph;  plastic  strain  was  derived  from 
these  values  by  subtracting  off  the  elastic  strain  as  indicated  by 
Hooke's  Law. 

Plastic  work  was  then  computed  from  a  graphical  presentation  of 
this  information.  A  trapezoidal  integration  was  done  on  the  graph  of 
stress  versus  plastic  strain  (see  Figure  12) .  The  value  of  plastic 
work  at  the  data  points  provided  the  abscissa  for  the  graph  representing 
the  relationship  of  Equation  (23). 


To  find  the  values  for  Z 


=  a  ■ 

i 


2D  ~\ 


2n  .  o 
n+1  ln73? 


,  the  value  of  plastic 


strain  rate  had  to  be  computed  for  each  data  point.  This  was  derived 
from  the  total  strain  rate  by  subtracting  off  the  elastic  strain  rate 
indicated  by  the  time  derivative  of  Hooke's  Law.  Then,  using  Equation 

m  i  * 

(24),  a  value  of  Z  was  computed  at  each  point.  The  data  was  presented 
graphically  as  in  Figures  13  through  16. 

One  problem  revealed  in  this  procedure  is  that  there  is  not  a  clear 
best  straight  line  approximation  to  the  data.  Perhaps  this  is  due  to 
scatter  of  the  experimental  data;  perhaps  the  model  simply  does  not 
adequately  represent  the  material  behavior.  Whatever  the  cause,  the 
data  used  in  this  study  seemed  to  follow  a  bi-linear  distribution,  and 


three  lines  were  taken  from  each  graph  (Figures  13  -  16)  as  representa¬ 
tive  of  the  data.  At  this  point,  for  each  of  the  four  values  of  n 


Figure  13.  Determination  of  m  and  Z 


Figure  16.  Determination  of  m  and  Zn,  n  =  3.211 


Some  parametric  studies  were  next  undertaken,  in  an  attempt  to 
refine  the  selection  of  values  for  the  coefficients  n  and  Z^.  A  set  of 
dimensionless  strain  rate  versus  stress  curves  were  prepared;  the  exper¬ 
imental  data  and  the  Bodner  model  predictions  were  plotted  for  each 
value  of  n  and  Z^.  The  results  are  contained  in  Figures  17  through  20. 
The  fit  of  model  to  data  in  each  case  is  measured  by  how  nearly  the  curve 
of  the  model  passes  through  the  two  (high  strain  rate,  no  recovery)  data 
points  indicated  by  arrows.  In  each  case,  the  fit  is  quite  good.  With¬ 
in  the  limits  on  accuracy  imposed  by  the  graphical  techniques  employed, 
all  values  are  arguably  adequate. 

Stress-strain  curves  were  used  next  to  compare  the  model  predictions 
with  the  data  in  another  fashion.  Values  for  A,  r,  and  Z^  were  selected 
without  too  much  concern  for  their  accuracy;  because  these  coefficients 
appear  only  in  the  recovery  term  of  Equation  (15)  ,  they  do  not  impact 
material  behavior  at  strain  rates  too  fast  for  recovery  to  be  signifi¬ 
cant.  This  is  illustrated  in  the  studies  of  Figures  48,  49,  51,  52,  54, 
and  55.  Using  the  different  values  of  n,  stress-strain  curves  were  then 
drawn  to  investigate  whether  varying  n  had  any  effeet  on  the  stress- 

strain  curve  of  the  model.  The  results,  contained  in  Figure  1,  show 

•  •  • 

that  the  four  curves  are  barely  distinguishable  from  one  another. 

Because  there  was  apparently  no  basis  in  the  data  for  discrimina¬ 
ting  between  the  four  values  of  n  and  Z^,  the  value  of  n  *  3.0  was 
selected  simply  for  the  convenience  provided  by  an  integer  value.  This 
choice  is  rather  arbitrary;  comparison  with  available  data  show  only 
that  the  model  predictions  based  on  n  =  3.0  fit  the  data  as  well  as  the 
predictions  based  on  any  other  n  value  under  consideration. 


NIVUIS 
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Figure  20.  Fit  of  Rodner  Model  to  IHpJi  Strain  Kate  Data,  n  -  3.211 


The  data  base  used  in  this  work  provided  only  two  points  (e  =  1.3  x 
-5  -1  •  -3  -1 

10  sec  and  e  =  1.6  x  10  sec  )  with  which  to  evaluate  n  and  Z^. 

•  -4  -1 

Having  a  third  data  point  —  at,  say,  £  =  1.0  x  10  sec  —  would  have 
allowed  significantly  more  confidence  in  the  values  thus  established  for 
n  and  Z^.  The  relatively  strain  rate  insensitive  behavior  of  IN  718 
under  the  conditions  of  these  tests  allowed  appropriate  coefficient 
values  to  be  determined  from  only  two  data  points.  For  a  more  strain 
rate  sensitive  material,  however,  the  exact  determination  of  coefficient 
values  would  likely  require  at  least  three  appropriate  data  points. 

With  a  value  of  n  =  3.0  and  the  corresponding  Z^  value  selected, 
the  material  behavior  predicted  by  the  Bodner  model  for  the  three  values 
of  m  (and  Z^)  were  compared  with  the  data.  It  is  apparent  from  these 
results,  presented  in  Figure  21,  that  the  choice  of  values  fcTr  m  and 
significantly  impact  the  shape  of  the  model's  predicted  stress-strain 
curve.  If  the  m-ZQ  values  suggested  by  the  steepest  line  fit  of  Figure 
14  (which  corresponds  to  m  =  .417  MPa  ^  and  Z^  =  1622  MPa)  are  used,  the 
front  of  the  model's  stress-strain  curve  in  Figure  21  fits  the  experi¬ 
mental  data  more  closely.  If  the  shallowest  line  fit  of  Figure  14  is 

used  (corresponding  to  m  *  .105  MPa  ^  and  Z^  =  1725  MPa),  then  the  back 

*  •  * 

of  the  stress-strain  curve  fits  the  data  more  closely.  And  if  the 
"mean"  line  of  Figure  14  is  used  (m  =  .219  MPa  ^  and  Z^  =  1669  MPa)  the 
resulting  stress-strain  curve  doesn't  quite  fit  the  data  either  at  the 
front  or  the  back,  but  it  isn't  too  far  off  in  either  place. 

The  values  of  m  and  Z^  which  enabled  the  model  to  closely  fit  the 
front  of  the  stress-strain  curve  were  selected  as  the  best  fit.  The 
material  behavior  in  the  smaller  strain  range  is  of  greater  engineering 
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interest;  in  cases  of  cyclic  loading  this  is  the  range  of  strains  that 
is  most  likely  to  be  studied.  Recognizing  that  material  behavior  pre¬ 
dicted  by  the  model  for  large  strains  will  not  likely  fit  experimental 
data  very  well,  values  of  m  =  .417  MPa-1  and  ZQ  =  1622  MPa  were  selected 
as  representative  of  the  material. 

After  values  were  determined  for  A,  r,  and  Z  the  matter  of  appro¬ 
priate  values  for  m  and  was  again  examined.  Analytical  predictions 
according  to  the  Bodner  model  of  material  behavior  in  creep  tests  were 
compared  with  experimental  data.  Coefficient  values  depending  on  the 
"front"  fit  of  Figure  14  (m  =  .417  MPa  1  and  Z^  =  1622  MPa)  and  values 
depending  on  the  "mean"  fit  (m  =  .219  MPa  1  and  Z^  =  1669  MPa)  were  both 
used.  The  results  of  this  study  are  depicted  in  Figures  22  and  23. 

Figure  22  shows  the  Bodner  model  predictions  under  a  135  KSI  load 
using  both  the  "front"  fit  values  of  m  and  ZQ  and  the  "mean"  fit  values. 
Under  this  load,  the  "mean"  fit  values  provide  a  very  good  model  of  the 
secondary  creep  strain  rate  —  better  than  that  provided  by  the  "front" 
fit  values  of  m  and  Z^.  Yet  Figure  23  —  which  depicts  the  same  infor¬ 
mation  as  in  Figure  22  except  the  load  is  140  KSI  —  shows  that  neither 
pair  of  values  gives  a  very  accurate  model  of  the  measured  material 
behavior.  The  inaccuracies  causfed  by  the  curve  fitting  techniques  on 
semi-log  axes  are  revealed  here;  an  error  in  predicted  strain  rate  of 
the  same  order  as  the  strain  rate  is  as  good  as  can  be  expected.  It 
would  appear  that  creep  behavior  does  not  lend  itself  to  use  in  evalu¬ 
ating  m  and  Z^  values. 

With  values  selected  for  n,  m,  Z^  and  Z^,  parametric  studies  were 
again  used  to  determine  appropriate  values  for  the  recovery  coefficients 


Figure  22.  Creep  Curves:  Comparison  of  Bodner  Model  With  Data 
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A,  r,  and  Z^.  The  general  behavior  of  A,  r,  and  Z 2  are  illustrated  in 
Figures  6,  7,  and  8.-  It  was  decided  to  fit  A  to  the  data  first,  then 
Z^,  and  lastly  r,  using  the  strain  rate  versus  stress  presentation  to 
establish  the  appropriate  coefficient  values. 

The  coefficient  A  controls  the  onset  of  significant  recovery 
effects.  Figure  24  presents  the  parametric  study  used  to  select  a  value 
of  A  =  . 0015  sec  ^ . 


The  coefficient  Z^  controls  where  the  very  low  strain  rate  material 
behavior  becomes  strain  rate  insensitive  once  more.  This  region  of  the 
strain  rate  versus  stress  curve  is  of  little  interest  in  practical 
engineering  applications,  and  without  a  significantly  more  extensive 
data  base  an  accurate  value  of  Z^  cannot  be  established.  For  the  pur¬ 
pose  of  these  applications  it  is  sufficient  that  the  return  CO  strain 
rate  insensitive  behavior  occur  below  the  strain  rate  range  of  the  data. 
Figure  25  contains  the  parametric  study  done  on  Z^  at  this  point;  a 
value  of  Z^  =  718  MPa  was  chosen.  If  a  smaller  value  had  been  selected, 
the  curve  would  have  indicated  a  return  to  strain  rate  insensitive  be¬ 
havior  at  an  even  lower  strain  rate.  This  would  not- have  affected  the 

ability  of  the  Bodner  model  to  predict  material  behavior.  The  only 

«  •  * 

effect  would  have  been  to  change  the  value  selected  for  r  —  since  Z ^ 
also  affects  the  slope  of  the  center  portion  of  the  curve. 

The  coefficient  r  controls  the  slope  of  the  mid-section  of  the 
strain  rate  versus  stress  curve.  Figure  26  contains  the  parametric 
study  used  to  arrive  at  the  value  r  =  7.0. 

At  this  point  a  complete  data  set  has  been  determined.  It  is 
tabulated  in  Table  IV.  The  fit  of  the  Bodner  model  to  the  experimental 


TABLE  IV 


Coefficient  Values  for  the  Bodner  Model,  IN  718,  1200°F 


Coefficient 

Description 

Value 

Do 

Limiting  value  shear  strain  rate 

in6  -1  i 

10  sec 

n 

Strain  rate  sensitivity  parameter 

3.0 

m 

Rate  of  work  hardening 

.417  MPa-1 

zo 

Initial  value  of  hardness 

1622  MPa 

Z1 

Saturation  hardness 

1795  MPa 

Z2 

Minimum  recoverable  hardness 

718  MPa 

r 

Exponent  on  recovery  term 

7.0 

A 

Coefficient  in  recovery  term 

.0015  sec  ^ 

E 

Young ' s  modulus 

162.5xl03  MPa 

Parametric  Study  Used  To  Evaluate 


.00  0.20  0.4  o  0-60  0-R0  1  -  00 

STRESS 

Figure  26.  Parametric  Study  Used  To  Evaluate 


data  using  these  values  is  presented  in  Figure  27. 

It  is  apparent  from  this  figure  that  r.o  attempt  was  made  to  bring 

the  curve  predicted  by  the  Bodner  model  through  the  two  lowest  data 

*  -9  -1  •  -8  -1 

points  (£  =  4.0  x  10  sec  and  e  =  1.6  x  10  sec  ).  The  model  lacks 

the  flexibility  required  to  bring  the  curve  back  up  through  those  two 
points  before  it  turns  downward  sharply  in  the  very  slow  strain  rate 
portion  of  the  curve.  The  very  slow  strain  rates  associated  with  these 
points  render  that  portion  of  the  data  less  important  to  engineering 
applications;  matching  the  rest  of  the  data  was  considered  adequate  at 
this  time. 

Thus  values  were  established  for  the  coefficients  of  the  Bodner 
model  for  Inconel  718.  Figure  28  serves  as  a  brief  descriptive  guide 
to  the  steps  followed  during  the  evaluation  process. 

Studies  of  material  behavior  under  cyclic  loading  conditions  were 
performed  using  the  coefficient  values  determined  in  this  work.  Various 
strain  rate  and  stress  control  situations  were  explored.  The  results  of 
this  study  are  discussed  in  the  following  sections. 


.125  .25  .375  .50  .625  .75  .875 

STRESS 

Figure  27.  Fit  of  the  Bodner  Model  to  Data,  All  Coefficients  Evaluated 


EVALUATE  N  AND  Z ,  USING  HIGH  STRAIN  RATE  DATA  AND  EQUATION  (21). 


REPEAT  EVALUATION  OF  A,  Z0  AND  R 


Load  Cyclins  Effects  Using  The 


Bodner  Constants  As  Evaluated 


Analytical  studies  were  conducted  of  the  material  behavior  under 
cyclic  loading  as  predicted  by  the  Bodner  material  model  developed  in 
previous  sections.  The  formulation  of  the  equations  incorporated  into 
this  study  is  contained  in  Equations  (16)  through  (20).  Fully  isotropic 
(q  =  1)  and  fully  directional  (q  =  0)  strain  hardening  conditions  were 
modeled.  Additionally,  Bodner  has  indicated  [7]  that  it  is  possible  to 
represent  strain  softening  by  incorporating  a  negative  value  of  q  into 
the  model.  This  last  case  is  interesting,  for  if  it  compares  well  with 
data  then  the  Bodner  material  model  is  very  comprehensive.  It  would 
show  the  model  can  consider  the  range  of  material  characteristics  from 
strain  hardening  to  softening,  including  anisotropic  plastic -hardening . 
The  results  of  the  uniaxial  analyses  were  compared  with  experimental 
fatigue  data.  This  section  describes  in  detail  the  cyclic  loading  study 
carried  out. 

Figures  29  and  30  contain  experimental  data  obtained  from  the  Air 
Force  Wright  Aeronautical  Laboratories  for  IN  718  at,  1200°F  gathered 
under  cyclic  loading  conditions.  Figure  29  depicts  the  material  be- 
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havior  measured  for  a  strain-controlled  test  (±  .8.%  strain)  while  Figure 
30  depicts  the  material  behavior  measured  during  a  stress-controlled 
test  (±  120  KSI) .  From  these  two  figures  it  is  clear  that  IN  718  ex¬ 
hibits  some  strain  softening  behavior.  Later  cycles  require  lower 
stress  levels  to  achieve  comparable  strains  than  do  earlier  cycles. 

This  phenomenon  is  especially  prominent  in  Figure  29,  where  the  stress 


level  for  maximum  strain  is  getting  smaller  in  absolute  value  for  each 


successive  cycle.  Figure  30  indicates  that  the  plastic  work  is  continu¬ 
ously  increasing  over  the  cycles  shown.  Both  of  these  phenomena  in  each 
figure  are  observed  over  a  limited  number  of  cycles  and  therefore  judge¬ 
ments  that  are  construed  are  based  entirely  on  this  limited  cyclic 
number.  Further  cyclic  loading  in  particular  may  reduce  the  amount  of 
strain  softening  evidenced  each  cycle.  Figures  31  and  32  illustrate 
the  behavior  predicted  by  the  Bodner  model  for  IN  718  using  values  of 
q  =  1  (Figure  31)  and  q  =  0  (Figure  32).  In  both  cases,  the  cycles 
repeat  without  change  as  early  as  the  third  cycle.  The  behavior  depicted 
in  these  figures  is  that  of  strain  hardening.  Figure  31  shows  isotropic 
hardening  (q  =  1)  while  Figure  32  shows  directional  strain  hardening 
(q  =  0).  Comparison  of  these  two  figures  with  the  experimental  data  of 
Figure  29  indicates  that  these  Bodner  model  predictions  represent  the 
first  cycle  fairly  well  but  do  not  indicate  the  softening  behavior 
evidenced  by  the  material.  This  was  to  be  expected;  values  of  q  between 
zero  and  one  imply  isotropic/kinematic  strain  hardening,  not  softening. 

It  was  suggested  [7]  that  q  <  0  be  used  to  model  cases  of  strain 
softening.  It  is  enlightening  to  examine  what  effect  a  negative  value 
—  say,  q  =  -.1  —  has  on  the  cyclic  equations  governing  changes  in 
nardness  in  this  formulation  (Equations  (16)  through  (20)). 

Initially,  as  Equation  (16)  indicates, 

ZC  =  ZC  =  Z0  (16) 

Given  that  the  initial  loading  is  tensile,  observe  that 

U  =  1  (28) 

a 
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Thus  the  two  time  rates  of  change  of  hardness  are 


ZC  =  -.1  Z  +  (1.1)  Z  = 

z<i) 

(29a) 

ZC  =  -.1  Z  -  (1.1)  Z  = 

-1-2  *u> 

(29b) 

where 

Z(1)  =  m  (ZL  -  Zt)0£P  -  A 

(30) 

At  the  end  of 

the  first  loading  time  increment  At,  the  two  values 

of 

hardness  are 

Z(l)  "  zo  +  z(i)Ati 

Ola) 

Z(l)  =  zo  "  1,2  Z(l)At3 

(31b) 

It  can  be  seen  from  this  that  while  the  hardness  of  the  specimen  in 

tension  is  increasing,  the  hardness  in  compression  is  decreasing  at  a 

faster  rate.  As  soon  as  the  load  becomes  compressive,  it  is  the  value 

of  Z  that  is  used  to  express  material  hardness.  Since  this  value  is 

less  than  the  value  for  Z ^  at  this  time  —  indeed,  less  than  the  original 

value  Z  —  it  is  clear  that  the  material  will  be  softer  in  compression 
U  .  s 

loading  subsequent  to  the  application  of  the  tensile  load. 

Figure  33  illustrates  this  behavior  through  a  graph  of  the  computer 
t  c 

output  of  Z  and  Z  over  eight  cycles  using  a  value  of  q  *  -.1.  The 
t  c 

values  of  Z  and  Z  were  plotted  every  quarter  cycle  to  show  the  general 

behavior  of  the  hardness  values  without  going  into  much  detail.  As  the 

value  of  Zt  increases  during  the  application  of  a  tensile  load,  the 
c  c 

value  of  Z  decreases;  as  the  value  of  Z  increases  during  the  application 


Behavior  of  Z  and  Z  During  Application 


of  a  compressive  load,  the  value  of  Z  decreases.  And  with  each  cycle, 

t  c 

the  average  value  of  Z  is  lower.  Each  time  Z  and  Z  reach  a  minimum 
value,  it  is  less  than  the  cycle  before;  each  time  they  reach  a  maximum 

value,  it  is  less  than  the  cycle  before.  This  decrease  in  Z  is  what  is 

called  strain  softening  and  it  can  be  directly  related  to  the  character 
of  the  stress-strain  curve. 

Figures  34  and  35  depict  the  Bodner  model  predictions  for  q  =  -.05 
and  q  =  -.1,  respectively.  As  in  the  cases  of  q  =  0  and  q  *  1,  the 

model's  loop  is  quite  similar  to  that  of  the  data,  both  in  size  and 

shape.  These  two  figures,  however,  possess  an  additional  characteris¬ 
tic;  they  show  clearly  the  strain  softening  trend  illustrated  by  the 
experimental  data.  The  loops  do  not  stabilize  in  these  early  cycles; 
the  material  softening  is  still  evidenced  by  both  model  and  data  in  the 
tenth  cycle.  The  similarity  between  Figure  29  and  Figures  34  and  35  is 
striking.  Small  negative  values  of  q  do  indeed  enable  the  modeling  of 
the  behavior  of  IN  718  under  a  cyclic  load. 

Consider  now  the  case  of  stress-controlled  test  situations.  Figure 
30  illustrates  the  experimentally  measured  behavior  of  IN  718  under  a 
cyclic  load  of  ±  120  KSI.  The  evidence  of  strain  softening  is  dramatic. 
Figures  36  and  37  contain  the  Bodner  model  predictions  for  q  *  1  and 
q  =  0;  not  surprisingly,  the  model  predictions  do  not  match  the  data 
very  well,  overall.  The  early  cycles  are  matched  fairly  well,  but  again 
the  model's  loop  stabilizes  where  the  experimental  one  continues  to  grow 

Figures  38  and  39,  however,  model  the  material  behavior  quite  well, 
using  values  of  q  *  -.05  and  q  =  -.1.  The  shape  of  the  initial  loading 
cycle  matches  the  data  in  both  size  and  shape.  The  model  then  goes  on 


Bodner  Model  Prediction  of  Cyclic  Loading  Behavior:  Strain  Control 


Bodner  Model  Prediction  of  Cyclic  Loading  Behavior:  Stress  Control 


Figure  38.  Bodner  Model  Prediction  of  Cyclic  Loading  Behavior:  Stress  Control, 


Bodner  Model  Prediction  of  Cyclic  Lending  Behnvior:  Stress  Co 


Co  predict  a  growing  and  spreading  of  the  loops  just  as  the  experimental 
data  shows.  The  loops  of  the  data  are  growing  faster  than  the  model's, 
but  the  trend  is  the  same. 

It  is  thus  seen  that  the  Bodner  material  model  possesses  a  remark¬ 
able  flexibility  in  the  realm  of  plasticity.  It  can  represent  both 
isotropic  and  directional  strain  hardening  as  well  as  strain  softening. 
With  the  addition  of  only  a  single  parameter,  q,  the  formulation  of 
Equations  (16)  through  (20)  allowed  the  extension  of  the  model  from 
monotonic  loading  to  cyclic  load  application,  with  remarkable  results. 

There  is  a  problem,  however,  that  limits  the  uses  of  this  model 
for  strain  softening  situations  with  negative  values  for  q.  There  is 
no  mechanism  in  this  formulation  to  impose  a  minimum  value  of  hardness. 
There  is  no  constraint  in  how  small  Z  may  get;  eventually  a  negative 
value  is  reached  and  the  equations  make  no  sense.  Figure  33,  which  is 

based  in  Bodner ' s  model  for  the  case  of  q  =  -.1,  shows  how  the  values 
t  c 

of  Z  and  Z  approach  zero  when  the  mean  curve  is  considered.  This 
basic  short-coning  may  be  overcome  by  the  addition  of  a  parameter  limit¬ 
ing  the  strain  softening  representation.  For  a  small  number  of  cycles, 
though,  and  a  q  whose  absolute  value  is  less  than  or  equal  to  .1,  the 

j  • 

model  can  be  used  to  illustrate  strain  softening  material  behavior. 

The  limitations  imposed  by  this  aspect  of  Bodner's  model  are  not 
great,  considering  the  remarkable  flexibility  it  still  has.  It  is  clear 
that  the  potential  is  here  for  a  tool  that  will  handle  a  wider  range  of 
strain  softening  behavior;  all  that  is  required  is  the  incorporation  of 
a  mechanism  for  imposing  a  limiting  value  of  softness. 


Conclus ions 


The  following  conclusions  are  made  regarding  the  nature  and  evalu¬ 
ation  of  the  Bodner  material  coefficients  and  the  use  of  this  formula¬ 
tion  to  model  cyclic  loading  conditions.  These  results  are  generaliza¬ 
tions  based  on  work  dene  with  IN  718  at  1200°F  and  on  a  review  of  the 
literature  concerning  modeling  of  rate-dependent  material  behavior: 

1.  Evaluation  of  the  Bodner  model  material  coefficients  requires 
the  following  minimum  experimental  data  base: 

a.  Two,  preferably  three,  data  points  for  the  strain-rate 
versus  stress  curve  where  test  conditions  are  high  strain  rate  with  no 
recovery. 

b.  At  least  one  stress-strain  curve  done  under  high  strain 
rate/no  recovery  conditions;  the  curve  should  extend  out  at  least  to 
where  the  plastic  strain  equals  the  elastic  strain. 

c.  Several  data  points  spread  across  the  lower  strain  rate 
portion  of  the  curve,  spanning  the  strain  rates  of  interest  to  the 
researcher . 

2.  Each  coefficient  of  the  Bodner  model  has  characteristic  effects 
on  material  behavior  .that  can  befdiscerned  individually. 

3.  A  systematic  procedure  can  be  established  for  determining 
appropriate  values  for  the  Bodner  model  coefficients;  it  is  contained 
in  this  thesis. 

4.  Evaluation  of  coefficient  values  is  in  some  significant  amount 

•  tive  in  nature.  Some  judgement  must  be  made  regarding  which  strain 
•prise  the  region  of  interest;  some  evaluation  of  what  consti- 
•  fit"  to  the  stress-strain  curve  (e.g.  in  the  evaluation  of 
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r.  and  Z  )  must  be  made. 

5.  Because  of  the  subjective  nature  of  coefficient  evaluation,  no 
two  researchers  are  likely  to  arrive  at  the  same  exact  coefficients 
from  a  given  data  base.  Differences  in  material  behavior,  if  any,  will 
be  traceable  to  the  fundamental  assumptions  of  the  individual  researcher 
(e.g.  which  strain  rates  are  important,  which  data  points  are  "bad  data", 
etc . )  . 

7.  Cyclic  behavior  predicted  under  this  formulation  of  the  consti¬ 
tutive  equations  models  experimental  results  quite  well.  Both  strain 
hardening  and  strain  softening  (for  sufficiently  small  values  of  q  and 
relatively  few  cycles')  can  be  modeled. 

8.  Further  development  of  the  cyclic  formulation  of  the  equations 
must  incorporate  a  mechanism  for  imposing  a  minimum  hardness'that  may  be 


reached  under  softening  conditions. 
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APPENDIX  A 


Discussion  of  the  Coefficients 

This  appendix  contains  a  detailed  discussion  of  each  material  co¬ 
efficient  of  the  Bodner-Partom  plastic  flow  law.  The  parametric  studies 
serve  to  illustrate  the  impact  of  the  various  coefficients  on  material 
behavior  predicted  by  the  Bodner  model.  The  figures  are  ordered  so 
that  the  variables  are  presented  in  the  order  in  which  they  were  studied 
in  the  applications  section.  Figure  5  provides  some  general  guidance 
regarding  the  information  contained  in  the  figures.  Each  figure  is 
accompanied  by  a  discussion  of  the  properties  revealed  by  the  graphical 
presentation. 

The  Coefficient  PQ 

Dq  is  the  limiting  value  of  the  plastic  strain  rate  in  shear.  For 

this  entire  exercise,  it  has  been  assumed  to  have  a  constant  value  of  10 

-1  4  . 

sec  .  Earlier  work  has  used  a  value  of  10  [1,2 J;  however,  recent  work 

[19]  indicates  that  a  value  of  10^  is  more  appropriate.  D^  effects  at 

«I.at  point  the  ultra  high  strain  rate  portion  of  the-strain  rate  versus 

stress  curve  bends  off  to  asymptotically  approach  D^;  for  as  strain 

m  •  * 

rates  approach  Dq,  the  stress  necessary  to  sustain-  these  strain  rates 
becomes  infinite.  Changing  the  value  of  D^  biases  the  entire  curve  up 
or  down,  but  does  not  effect  the  relationships  between  the  other  co¬ 
efficients. 

The  Coefficients  n  and 

The  coefficient  n  describes  the  material's  strain  rate  sensitivity. 
The  high  strain  rate  version  of  Equation  (10)  (Equation  (21))  contains 


Thus  n 


only  —  assumed  to  be  constant  at  10  sec  —  n,  and  Z^. 
and  control  a  significant  portion  of  the  material's  high  strain  rate 
behavior.  The  constant  n  has  been  determined  [18]  to  be  inversely 
related  to  temperature;  the  form  of  the  relationship  is 

n  =  f  +  b  (32) 

The  coefficient  n  is  also  inversely  related  to  the  "intrinsic  velocity 
of  dislocation  motion"  [l]  in  the  atomic  lattice  —  meaning  that  lower 
values  of  n  correspond  to  higher  dislocation  velocities  and  freer  move¬ 
ment  of  dislocations.  The  coefficient  n  is  not  dependent  on  loading 
history  [l].  It  is  thus  descriptive  of  the  way  the  material  itself 
responds  to  the  introduction,  propagation,  and  elimination  of  disloca¬ 
tions,  not  to  the  effect  those  dislocations  have  on  material  "behavior. 

Values  for  n  and  Z^  are  tied  together  mathematically;  at  high  strain 
rates  the  relationship  is  that  of  Equation  (21): 


/OpP  , -i 

In  (-In  (-—-))  =  -2n  In  a  +  (2n  In  Z.  +  In  (~±)) 
2D0  1  2n 


Z^  represents  the  saturation  (maximum)  value  of  Z.  In  a  high  strain 
rate  situation,  wher£;  there  is  no  recovery,  Z  =  Z^. 


Figure  40 

This  figure  depicts  the  combined  effects  of  variations  in  n  and  Z^. 
Strain  rate  sensitivity,  dependent  on  n,  is  observed  in  that  high  (but 


not  ultra  high)  strain  rate  section  of  the  curve  where  recovery  effects 


are  negligible.  Decreasing  n  has  the  effect  of  increasing  strain  rate 
sensitivity;  this  is  seen  in  the  shallower  slope  of  the  curve  at  lower 
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Figure  40.  Variation  of  n,  All  Other  Coefficients  Constant 


values  of  n.  As  the  strain  rates  go  down,  recovery  effects  are  no 

m  • 

longer  negligible,  and  the  significance  of  a  given-value  of  n  as  a  strain 
rate  sensitivity  parameter  becomes  obscured.  At  high  strain  rates, 
where  Z  =  Z^,  the  equations  reduce  to  the  closed  form  of  Equation  (24): 

Z.  /  o  f  2Dn  lYsn 

—  ~(-zrr  In  — -  )  (24) 

0  V"+1  I-  /3£PJ/ 

Again,  D  (a  constant),  n  and  Z1  are  the  only  coefficients  present;  hence 


n  and  are  Che  only  coefficients  capable  of  modifying  the  high  strain 
rate  portion  of  the  strain  rate  versus  stress  curves  of  this  parametric 
study. 

The  translation  of  the  curve  along  the  x-axis  is  due  to  the  choice 
of  as  the  non-dimensionalizing  agent. 

Figures  41  and  42 

The  range  of  values  of  n  and  which  emerge  from  the  interpretation 

of  the  data  (Figures  9  and  56)  is  presented  in  Figures  41  and  42.  In 

Figure  41  there  is  clearly  more  than  one  saturation  stress  represented. 

This  is  due  to  the  fact  that  the  point  B  of  Figure  56  (through  which  the 

four  lines  representing  the  values  for  n  all  pass)  corresponds  to  a 

-3  -1 

strain  rate  of  1.185  x  10  sec  .  Only  at  that  strain  rate  do  all  four 
n  values  correspond  to  the  same  saturation  stress.  Yet  the  variation  in 
n  from  3.0  to  3.211  is  not  significant  in  Figure  41.  From  Figure  42  it 
is  clear  that  the  change  in  the  curve  slope  caused  by  changing  n  such  a 
small  amount  is  not  significant.  An  arbitrary  choice  of  n  =  3.0  was 
made,  based  on  a  graphical  comparison  with  experimental  data  (Figures 
16  through  19) . 

•  t 

The  Coefficients  m  and  Z^ 

The  coefficient  m  describes  the  rate  of  work  hardening  exhibited  by 
the  material.  Its  value  also  determines  how  long  work  hardening  effects 
will  be  dominant  and  recovery  effects  negligible  as  strain  rate  is 
decreased.  Mathematically,  it  is  the  value  of  m  coupled  to  the  value  of 
A  that  determines  above  what  strain  rates  hardening  effects  dominate  so 


completely. 


Variation  of 


Zq  is  linked  mathematically  to  m  much  as  is  linked  to  n.  This 
relationship  is  expressed  in  Equation  (23): 


In  (Z1  -  Z)  =  In  (Z1  -  ZQ)  -  mWp  (23) 

Zq  represents  the  initial  state  of  hardness  in  the  present  loading;  the 
lower  the  value  of  Z^,  the  more  hardening  must  be  accomplished  before 
saturation  is  reached  and  Z  =  Z^.  Z^  is  dependent  on  temperature. 


Figure  43 

The  coefficient  m  times  Z^  is  the  non-dimensional  form  of  the  co¬ 
efficient  describing  the  rate  of  work  hardening,  m.  Once  a  value  is 
established  for  Z^,  the  behavior  depicted  by  Figure  43  represents  vari¬ 
ations  in  m.  If  one  compares  Figure  43  with  Figure  49,  considering 
Equation  (15): 


(15) 


it  is  observed  that  A  and  m  hold  analagous  positions  as  multiplicative 
coefficients  in  the  two  terms  of  Equation  (15).  Increasing  A  has  much 
the  same  effect  as  decreasing  m,  which  is  to  evidence  the  effects  of 
recovery  phenomena  at1  higher  stra’in  rates.  The  interface  between  the 
high  strain  rate  regime  (where  there  are  no  significant  recovery  effects) 


and  the  strain  rate  region  where  both  recovery  and  hardening  are 
important  is  controlled  by  this  balance  between  A  and  m. 


Variation  of  m.  All  Other  Coefficient 


Figure  45 

The  coefficients  m  and  Z^  are  linked  mathematically.  Lower  values 
of  m  are  paired  with  higher  values  of  Z  .  Figure  45  depicts  the  effect 
of  choosing  various  interpretations  of  the  data  graphed  in  Figure  13. 
Using  the  m  and  values  from  the  first  linear  portion  of  the  data  of 
Figure  13  (the  "front"  values  of  Z^  and  m)  results  in  a  better  fit  of 
the  front  of  the  stress-strain  curve  of  the  data.  Using  the  m  and 
values  from  the  second  linear  portion  of  the  data  of  Figure  13  (the 
"back"  values)  results  in  a  flatter  stress-strain  curve  that  matches  the 
saturation  stress  level  of  the  data  more  precisely.  The  choice  of  m  and 
ZQ  significantly  effects  the  ability  of  the  model  to  accurately  predict 
material  behavior  as  depicted  on  the  stress-strain  curve. 

Figure  46 

The  effect  of  varying  is  illustrated  in  Figure  46.  Since  Zq 
represents  the  state  of  hardness  at  the  beginning  of  the  present  loading, 
the  difference  between  the  hardness  at  the  start  and  the  maximum  hardness 
Z^  describes  the  hardening  that  may  go  on  during  the  loading. 

The  Coefficient  A 

#  '  * 

A  is  the  coefficient  in  the  expression  for  the  thermal  recovery  of 
work  hardening.  The  magnitude  of  A  (relative  to  the  magnitude  of  m) 
determines  at  which  strain  rates  recovery  becomes  significant.  For  lower 
values  of  A,  thermal  recovery  remains  a  less  significant  phenomenon  until 
relatively  slower  strain  rates  are  reached.  Changes  in  A  do  not  signifi¬ 
cantly  impact  material  behavior  at  strain  rates  too  fast  for  recovery 
effects  to  be  significant,  since  A  is  tied  strictly  to  the  recovery  mech¬ 
anism. 
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Figure  47.  Variation  of  A/DQ,  All  Other  Coefficients  Constant 
Figure  47 

a  >  •  ' 

The  coefficient  A  (or  A/Dq  in  the  non-dimensionalized  form)  is  the 
coefficient  in  the  expression  for  thermal  recovery  of  work  hardening. 

The  "branching"  effect  in  the  lower  left  of  this  figure  shows  how  the 
primary  effect  of  changing  A  is  to  extend  the  region  where  recovery  is  a 
significant  feature  of  material  behavior.  The  onset  of  thermal  recovery 
is  determined  by  this  variable  (once  a  value  for  m  is  established  througt 
study  of  the  high  strain  rate  material  behavior);  for  larger  values  of  A 


thermal  recovery  dominates  at  higher  strain  rates. 


Figures  48  and  49 

Changes  in  A  are  not  significant  for  strain  rates  faster  than  that 

at  which  recovery  becomes  evident.  Figure  48  is  an  analytically  gener- 

•  -5  -1 

ated  stress-strain  curve  for  e  =  1.3  x  10  sec  ;  it  shows  clearly  that 

changes  in  A  for  A  at  least  as  small  as  0.001  are  not  significant  for 

this  strain  rate.  Figure  49  is  an  analytically  generated  stress-strain 

-6  -1 

curve  for  £  =  1.1  x  10  sec  ;  here,  changes  in  A  are  significant  until 
A  is  at  least  as  small  as  0.00001.  These  figures  show  that  it  is  a 
balance  between  the  strain  rate  in  question  and  the  magnitude  of  A  that 
determines  whether  smaller  values  of  A  will  impact  material  behavior  at 
the  strain  rate  in  question.  Larger  values  of  A,  whose  effect  is  to 
extend  the  region  where  recovery  effects  are  significant,  will  cause 
visible  changes  in  the  curve  as  soo«  as  the  implied  change  in  A  makes 
the  strain  rate  relatively  slow  enough  to  allow  recovery  effects  to 
surface. 


The  Coefficient  Z„  _ 

The  coefficient  Z ^  represents  the  material  hardness  at  the  completely 

g  ;  • 

non-work  hardened  —  or  fully  recovered  —  condition.  It  is  the  minimum 
recoverable  hardness  of  the  material,  and  it  is  dependent  on  the  temper¬ 
ature  conditions. 


Figure  50 

Once  a  value  is  established  for  Z^  through  study  of  the  material's 
high  strain  rate  behavior,  the  material  behavior  depicted  in  Figure  50 


'  becomes  that  of  Z.. 
*  2 
* 

> 


A  higher  value  of  Z7  represents  a  harder  condition 


CC'OOl 


TRA  I 


'  r\ 


l 


.10  .30  .50  .70  .90  1.10  130 

Figure  50.  Variation  of 

when  minimum  hardness  is  achieved  under  recovery.  On  Figure  50, 

•  i  • 

controls  where  the  curve  turns  downward  again  at  very  low  strain  rates. 
Figures  51  and  52 

Like  A,  Z ^  is  part  of  the  equation  descriptive  of  recovery  pheno¬ 
mena.  Figure  51  shows  that  changing  Z^  values  at  a  high  strain  rate, 
where  recovery  is  negligible,  has  no  significant  effect  on  the  material 
behavior  depicted  by  the  stress-strain  curve.  A  stress-strain  curve  set 
drawn  for  a  slower  strain  rate,  where  recovery  effects  become  significant 


is  shown  in  Figure  52.  Here  changes  in  are  clearly  visible. 

The  Coefficient  r 

The  coefficient  r  appears  as  an  exponent  in  the  expression  for 
thermal  recovery  of  work  hardening.  If  A  (or  A/Dq)  were  to  be  written 
as  a  function  of  the  onset  of  recovery,  r  would  be  the  curvature  of  that 
function.  Where  the  value  of  A  determines  at  which  strain  rate  recovery 
begins  to  be  important,  r  describes  how  quickly  recovery  rises  to  domi¬ 
nance.  As  with  A  and  Z the  effect  of  changing  r  is  not  seen  at  strain 
rates  too  fast  to  allow  thermal  recovery  of  plastic  work  to  be  signifi¬ 
cant. 

Figures  53,  54,  and  55 

Figure  53  shows  that  the  effect  of  changing  r  is  to  change  the 
sensitivity  of  the  recovery  mechanism  to  strain  rate.  Figures  54  and  55 
show  quite  clearly  that  changes  in  r  do  not  effect  material  behavior  at 
strain  rates  where  thermal  recovery  is  negligible.  At  a  strain  rate 
fast  enough  for  recovery  to  be  negligible  (Figure  54)  ,  changes  in  r 
are  not  distinguishable  on  a  stress-strain  curve.  At  a  slower  strain 
rate  where  recovery  effects  become  significant  (Figure  55),  increases 
in  r  are  reflected  iosa  higher  stress  level  of  the  stress-strain  curve 


over  the  domain  shown. 


Variation  of  r  at  a  Low  Strain  Rate 


APPENDIX  B 


Determination  of  n  and  Z^ 


When  both  stress  and  strain  rate  are  constant,  the  governing  equa¬ 
tion 


(10) 


says  that  Z  must  also  be  constant.  Further,  for  high  strain  rate  tests 
with  no  recovery  the  material  must  be  in  a  fully  work  hardened  state, 
"to  obtain  the  maximum  value  of  stress"  [9],  so  that  Z  =  Z^.  For  these 
steady  flow  conditions.  Equation  (21)  follows: 


In  -In 


A? 

2D„ 


=  -2n  In  O  +  (2n  In  +  In 


n+1 

2n 


(21) 


Thus,  it  is  seen  that  a  linear  relationship  exists  between  ln(-ln 


/3ep 

2D. 


) 


and  In  a.  If  high  strain  rate  experimental  data  is  graphed  on  appropri- 

/3ep 


ate  axes  (ln(-ln 


2D, 


)  versus  In  a)  the  slope  of  a  linear  approxima¬ 


tion  to  the  data  will  determine  a  value  for  n.  A  value  for  can  then 
be  computed  mathematically. 

The  test  data  used  for  this  ^procedure  consisted  of  three  constant 
strain  rate  tensile  tests  presented  in  Table  V: 


TABLE  V 


Experimental  Data,  Constant  Strain  Rate  Tensile  Tests 


Data  Point 

Strain  Rate 

Stress  (KSI) 

Stress  (MPa) 

1 

1.1x10  ^sec  ^ 

134 

923.906 

2 

1.3x10  ^sec  ^ 

142 

979.065 

3 

1.6x10  "^sec  ^ 

147 

1013.539 
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Data  point  1,  with  a  strain  rate  of  1.1  x  10  sec  ,  was  considered  at 
this  stage  even  though  the  strain  rate  is  suspected  of  being  slow  enough 
to  allow  non-negligible  recovery  effects  to  influence  the  data.  Figure 
56  shows  the  location  of  this  data  point  on  strain  rate  versus  stress 
axes.  If  straight  lines  are  drawn  to  characterize  the  two  upper  por¬ 
tions  of  the  data  (as  in  Figure  5),  data  point  1  is  found  to  lie  in  the 
strain  rate  region  where  recovery  effects  are  significant. 

Using  Equation  (21)  to  compute  the  values  for  graphing  leads  to  the 
following: 

TABLE  VI 

Experimental  Data,  Formatted  to  Allow  Determination  of  n  and 


Data  Point 

/3ep 

y  -  ln(-ln(  2J)  )) 

x  =  In  a 

1 

3.3207 

6.8286 

2 

3.2376 

6.8866 

3 

3.0154 

6.9212 

This  data  was  graphed  in  Figure  57.  It  was  clear  that  all  three  points 
were  not  part  of  a  single  linear  trend.  The  validity  of  point  1  was 
suspect  from  the  beginning;  that  point  was  simply  discarded.  The  strain 
rate  of  1.1  x  10  ^  sec  ^  was  apparently  too  slow  to  allow  recovery 
effects  to  be  neglected.  This  was  verified  by  parametric  studies  con¬ 
tained  in  Appendix  A,  after  a  complete  set  of  coefficients  was  obtained 
(Figures  48,  49,  51,  52,  54,  and  55).  At  a  strain  rate  of  1.1  x  10  ^ 
sec  \  effects  of  varying  A,  r,  and  Z ^  are  clearly  seen;  at  a  strain 
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Figure  56.  Experimental  Data  Base;  IN  718,  1200°F 
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POINT  B 
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rate  of  1.3  x  10  sec  they  are  hardly  clscernable.  Since  these  three 

coefficients  appear  only  in  the  recovery  term  of  Equation  (15)  these 

results  indicate  that  recovery  effects  are  not  negligible  at  e  =  1.1  x  10 
-1 

sec 

Working  with  the  remaining  two  data  points,  the  value  of  n  was 
computed  in  a  variety  of  ways.  First,  working  mathematically,  a  value 
of  n  =  3.211  was  computed  using  the  coordinates  of  the  two  points.  In 
order  to  illustrate  the  limitations  on  accuracy  imposed  by  the  graphical 
technique  that  would  have  been  used  had  more  than  two  appropriate  data 
points  been  available,  a  second  value  was  computed.  Coordinates  of  two 
points  (other  than  the  data  points)  were  picked  off  a  line  drawn  through 
the  data,  and  they  indicated  a  value  of  n  =  3.125.  Because  n  *  3.0  is  a 
nice  integer  value,  the  line  representing  it  was  also  drawn;  “it  appears 
to  also  be  fairly  representative  of  the  data.  As  a  lower  limit,  the 
line  representing  n  =  2.5  was  also  drawn.  These  four  values  represent 
the  range  of  n  values  to  be  considered. 

Values  for  Z^  were  then  calculated.  All  of  the  aforementioned  lines 
were  drawn  through  a  single  point  B  (whose  coordinates  are  (6.92,  3.03)) 
in  Figure  57.  Using  this  point,  the  various  values  for  n  and  Equation 

mi  • 

(16)  lead  to  the  following  values  for  Z^: 

TABLE  VII 

V  ues  Computed  for  n  and  Z, 


n 

Z1(MPa) 

3.211 

i732.871 

3.12.5 

1756.858 

3.0 

1794.672 

2.5 

1992.878 

Z1  (KSI) 


251.329 

254.808 

260.293 

289.040 


One  consequence  of  the  use  of  point  B  should  be  mentioned.  It  is 

/3rP 

only  at  the  strain  rate  indicated  by  ln(-ln(— -) )  =  3.03  that  the 
saturation  stresses  will  be  the  same  for  all  the  values  of  n.  It  is 
partly  because  of  this  fact  that  the  saturation  stress  are  not  the  same 
in  the  stress-strain  curves  where  n  is  varied  (e.g.  Figure  41). 


APPENDIX  C 


In  (Z1  -  Z)  -  In  (Zx  -  ZQ)  -  mWp  (23) 

Appropriate  high  strain  rate  data,  graphed  on  In  (Z  -  Z)  versus  W  axes 

l  P 

may  be  linearly  approximated  to  determine  m  from  the  slope  of  the  line 
and  Zg  from  the  y- intercept. 

The  data  used  for  this  procedure  came  from  a  single  high  strain 
rate  tensile  test  and  is  presented  in  Figure  58.  This  data  was  chosen 
because  the  strain  rate  was  high  enough  that  recovery  effects  would  most 
likely  be  negligible,  and  the  data  maintained  a  nice,  smooth  character 
across  a  wide  range  of  strains. 

The  elastic  modulus  [10]  was  first  drawn  on  the  figure.  It  happens 
that,  although  the  data. is  reasonably  accurate  in  the  inelastic  region, 
the  accuracy  in  the  elastic  region  is  not  good.  Using  the  elastic 
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modulus  as  a  more  reliable  description  of  the  material  s  elastic  be¬ 
havior,  the  modulus  was  drawn  on  the  figure,  tangential  to  the  curve  of 
the  data.  The  origin  of  the  axes  was  then  translated  a  negative  0.000575 
units  along  the  x-axis.  This  allowed  the  curve  to  indicate  zero  strain 
at  a  condition  of  zero  stress. 

Coordinates  of  selected  data  points  were  then  read  off  the  figure. 
The  following  table  contains  the  information  as  digitized  at  this  point: 


TABLE  VIII 

Experimental  Data  Frop  Figure  58 


Point  No. 

Stress (y) 

Strain (x) 

Strain  plus 
"offset"  strain 

0 

122.6 

KSI 

.00465 

.005225 

1 

126 

KSI 

.00492 

.005495 

2 

128 

KSI 

.00506 

.005635 

3 

130 

KSI 

.00535 

.005925 

4 

132 

KSI 

.00587 

.006445 

5 

134 

KSI 

.00654 

.007115 

6 

136 

KSI 

.00725 

.007825 

7 

433 

KSI 

.  .0085 

.009075 

8 

138.6 

KSI 

.00955 

.010125 

9 

139 

KSI 

.01114 

.011715 

10 

139.5 

KSI 

.012985 

.01356 

Computation  of  the  plastic  work  W  requires  that  plastic  strains  be 

P 

derived  from  the  total  strains  indicated  from  Figure  58.  Using  Hooke's 
Law,  plastic  strains  may  be  computed  from  the  equation 


The  resulting  values  were  then  plotted  against  stress  a,  and  a  numerical 
integration  technique  employed  to  evaluate  the  area  under  the  curve. 


TABLE  IX 

Computation  of  Plastic  Strain 


Point  No. 

Total  strain 

Elastic  strain 

Plastic  strain 

0 

.005225 

122. 6/23. 6xl03 

3.01  x  10“5 

1 

.005495 

126/23. 6xl03 

1.56  x  10  4 

2 

.005635 

128/23. 6xl03 

2.11  x  10~4 

3 

.005925 

130/23. 6xl03 

4.17  x  10-4 

4 

.006445 

132/23. 6xl03 

8.52  x"10'4 

5 

.007115 

134/23.6x10 3 

1.44  x  10-3 

6 

.007825 

136/23. 6xl03 

2.06  x  10-3 

7 

.009075 

138/23 . 6xl03 

3.23  x  10~3 

8 

.010125 

138. 6/23. 6xl03 

4.25  x  10“3 

9 

.011715 

139/23. 6xl03 

-  5.83  x  10"3 

10 

.01356 

139.5/23. 6xl03 

7.65  x  10~3 

i  • 


The  plastic  strain  indicated  for  the  point  //0  was  supposed  to  be 
zero,  and  is  included  here  to  indicate  the  limits  on  accuracy  imposed 
by  the  method  employed.  Linearly  interpolating  between  the  points  of 
Figure  12,  and  extrapolating  to  find  a  point  where  plastic  strain  is 
zero,  the  values  for  plastic  work  were  then  computed.  The  general 
form  of  the  equation  used  is 


^  .  ”  -A  ^  A  ™  ™  . 


*  v"-1  4  a""1  (£P(n)-  eP^'1))  +I(eP(n>-  £P(n_1))  (an_  {J(n-l)) 

p  2  ' 

(34) 

0  +  121. 79(. 0000301  -  0)  4  ^00000301  -  0)(122.6  -  121.79) 

0  +  .00367  +  Js(. 0000301)  (122. 6  -  121.79) 

.00368 

W°  +  122. 6(. 000156  -  .0000301)  +  %(. 000156  -  .0000301) (126-122 

.00368  +  .0154  +  .000214 

.01933 

wj  +  126(. 000211  -  .000156)  +  Jf(. 000211  -  . 000156) (128-126) 

.01933  +  .00696  +  .0000553 

.02635 

+  128 (.000417  -  .000211)  4  hi- 000417  -  .000211) (130-128) 
.02635  +  .026327  4  .000205 
.05282 

4  130(. 000852  -  .000417)  4  hi- 000852  -  . 000417) (132-130) 
.05282  4  .05658  4  .000435 
.1098 

4  132 (,,Q0144  -  .000852)  4  hi- 001434  -  . 000852)  (134-132) 
.1098  4  .07725  4  .000586 
.1877 

W5  4  134(. 00206  -  .00144)  4  4(. 00206  -  . 00144) (136-134) 

P 

.1877  4  .0838  4  .000625 
.2721 
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W  =  W+  136(. 00323  -  .00206)  +  4(. 00323  -  . 00206) (138-136) 

P  P  ' 

=  .2721  +  .1585  +  .00117 
=  .4317 

W8  =  W7  +  138(. 00425  -  .00323)  +  %(. 00425  -  . 00323) (138. 6-138) 

P  P 

*  .4317  +  .1414  +  .000307 
=  .5734 

W9  =  W9  +  138. 6(. 00583  -  .00425)  +  %(. 00583  -  .00425) (139-138.6) 
P  P 

=  .5734  +  .218  +  .000315 
=  .7918 

w10-  w9  +  139(. 00765  -  .00583)  +  *s(. 00765  -  .00583X139.5-139) 

P  P 

■=  .7918  +  .2535  +  .000456 
=  1.0457 


Thus  the  x-coordinates  of  Equation  (23)  were  computed.  Determina¬ 
tion  of  the  y-rcoordinates  begins  with  the  evaluation  of  the  equation 


Z 
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which  requires  that  plastic  strain  rates  be  determined  at  each  data 

«  » • 

point.  The  intervals  between  the  points  were  numbered  (0)  through  (10) 
Since  numerical  methods  were  used  to  determine  elastic  stress  rates 
(from  changes  in  the  stress  over  time  intervals) ,  strain  rates  for  the 
intervals  were  computed  and  the  values  of  two  adjoining  intervals 
averaged  to  arrive  at  an  appropriate  value  for  the  points. 

Plastic  strain  rates  were  determined  according  to  the  following 


relationships: 


e  *  1.3  x  10  *  c  +  e 


•e  _  o  _  Ao  1 
Z  “  E  At  E 


Since  the  total  strain  rate  is  constant,  t  can  be  computed  for  each 
interval  using  the  change  is  total  strain. 


TABLE  X 


Determination  of-  the  Elastic  Strain  Rate  e 


Interval 


0 


.000030085 

.000125932 

.000055254 

.000205254 

.000435255 

.000585254 

.000625254 

.001165254 

.00102447J 

.001573050 

.001823814 


At  (sec) 

Aa  (KSI) 

Aa  KSI 

At  SEC 

ec  (sec  1) 

2.314 

.8123 

.351 

1.49xl0“5 

9.687 

3.4 

.351 

1.49xl0~5 

4.2503 

2.0 

.471 

1. 99xl0~5 

15.789 

2.0 

.127 

5. 37xl0-6 

33.481 

2.0 

.0597 

2.53xl0-6 

45.019 

2.0 

.0444 

1.88xl0~6 

48.096 

2.0 

.0416 

1.7 6x1 0~6 

89.635 

2.0 

.0223 

9. 45x1 0"7 

78.814. 

0.6 

.00761 

3.23xl0~7 

121.004 

0.4 

.00331 

1.40xl0~7 

140.293 

0.5 

.00356 

1.51x1 0-8 

Values  for  the  plastic  strain  rates  over  the  intervals  can  be 
easily  computed,  as  the  total  strain  rate  is  constant  and  the  elastic 
component  has  now  been  calculated.  The  plastic  strain  rates  at  the 
points  were  then  determined  by  averaging  the  rates  over  the  intervals 


Interval 


0 


•  P  — 

er  (sec  )  Point 


-1.87x10 


E13  (sec  1) 

Stress (KSI) 

Stress(MPa) 

-1.87x10“ 6 

122.6 

845.305 

-4.4lxl0-6 

126 

868.748 

•  3.47x10“' 

128 

882.537 

9 . 05xl0-6 

130 

896.327 

1.08xl0“5 

132 

910.117 

1.12xl0“5 

134 

923.906 

1 . 16xl0-5 

136 

937.696 

1.24x10“ 5 

138 

951.933 

1.28xl0“5 

138.6 

954.933 

1.29xl0“5 

139 

958.380 

1. 28xl0-5 

139.5 

961.828 

Notice  that  the  plastic  strain  rates  over  intervals  (0),  (1),  and 
(2)  turned  out  to  hav,e  negative  values.  It  would  appear  that  the  points 
defining  these  intervals  are  sufficiently  close  to  the  elastic  region 
that  graphical  methods  are  not  accurate  enough  to  allow  adequate  deter¬ 
mination  of  the  strains  and  stresses  there.  Because  of  the  negative 
values,  the  data  corresponding  to  points  0  and  1  cannot  be  graphed,  and 
the  data  corresponding  to  point  2  is  highly  suspect  (and  will  be  ignored 
when  the  data  is  graphed)  because  one  of  the  neighboring  intervals  had  a 
negative  plastic  strain  rate. 


ffl 


TABLE  XIII 


TABLE  XV 


Computation  of  Ln(Z^-Z)  ;  n  =  3.211,  Z^  =  1732.871  MPa 


Z  (MPa) 

In  ( Z 1  -  Z) 

1.8025 

1590.744 

4.957 

1.7691 

1585.683 

4.992 

1.7672 

1608.346 

4.825 

1.7668 

1632.365 

4.610 

1.7664 

1656.310 

4.338 

1.7657 

1680.048 

3.967 

1.7654 

1685.802 

3.852 

1.7653 

1691.817 

3.715 

1.7653 

1697.908 

3.554 

This  data  is  presented  in  graphical  form  in  Figures  13  through  16. 
According  to  Equation  (23),  the  data  should  have  indicated  a  straight 
line  that  would  allow  determination  of  m  and  Z^.  In  each  case,  however 
the  data  appeared  to  follow  more  of  a  bi-linear  distribution  than  a 
linear  one.  When  the  data  point  2  (which  contains  some  influence  from 
that  negative  strain  rate)  is  ignored,  the  first  five  points  suggest 
one  linear  approximation  and  the  last  four  points  suggest  another.  It 
is  perhaps  a  curious  coincidence  that  the  distribution  is  bi-linear; 
the  failure  of  the  data  to  be  linear  could  be  due  to  the  inadequacy  of 
the  Bodner  model  or  simply  due  to  scatter  of  the  experimental  data. 
Whatever  the  cause,  one  single  straight  line  was  not  indicated  by  the 
data,  so  a  slight  modification  to  the  procedure  was  indicated. 

In  each  case,  three  straight  lines  were  taken  from  the  figure.  On 
straight  line,  drawn  through  the  first  five  data  points,  indicated  what 
is  referred  to  as  the  "front"  values  of  m  and  Z^.  One  straight  line, 
drawn  through  the  last  four  data  points,  indicated  what  is  referred  to 
as  the  "back"  values  of  m  and  Z^.  And  an  approximation  of  a  straight 
line  through  all  the  data  provided  what  is  called  tTTe  "mean"  values  of 
m  and  Z^. 

«  •  • 

The  following  are  the  four  data  sets  at  the  conclusion  of  this 
procedure: 

Set  #1:  n  **  2.5 

Z  =  1992.878  MPa 


"front" 

"back" 

mean 

m  =  .483  MPa 

-1 

m  =  .15  MPa  ^ 

m  «  .242  MPa 

-1 

Z„  =  1811.606 

MPa 

Z.  =  1924.161  MPa 

Z-  «=  1866.408 

MPa 

Set  “2 : 

n  =  3.0 

Z  =  1794.672  MPa 

"front" 

"back" 

mean 

m  =  .417  MPa  1 

m  =  .104  MPa-1 

m  *  .219  MPa-1 

ZQ  =  1622.24  MPa 

ZQ  =  1725.264  MPa 

ZQ  *  1669.461  MPa 

Set  03: 

n  *=  3.125 

Z  =  1756.858  MPa 

"front" 

"back" 

"mean" 

m  =  .4  MPa  1 

m  =  .0985  MPa-1 

m  *  .18  MPa-1 

ZQ  =  1586.143  MPa 

ZQ  *  1686.753  MPa 

ZQ  -  1637.754  MPa 

Set  "4: 

n  ■  3.211 

Z1  *  1732.871  MPa 

"front" 

"back" 

"mean" 

m  *  .4  MPa  1 

m  ■  .0938  MPa  1 

m  -  .185  MPa-1 

Zn  =  1562.155  MPa 

Zn  =  1663.463  MPa 

Z.  -  1608.284  MPa 

APPENDIX  D 


Computer  Programs 

This  appendix  contains  the  three  computer  codes  used  in  the  course 
of  this  work.  All  three  use  the  Bodner  constitutive  relations  to  pre¬ 
dict  material  behavior  for  uniaxial  loading  conditions.  The  two  programs 
called  BODNER  and  CYBODR  were  written  by  A.  M.  Rajendran  of  the  Univer¬ 
sity  of  Dayton  Research  Institute. 

Program  BODNER 

This  program  is  used  to  generate  stress- strain  curves  (e.g.  Figure 
1)  and  creep  curves  (e.g.  Figure  22).  It  is  written  in  FORTRAN  4;  it 
utilizes  an  external  data  file  (TAPE5)  for  input  and  creates  two  output 
files  (TAPE6  and  TAPE7).  The  file  TAPE6  contains  the  incremental  values 
of  the  program  parameters  and  is  formatted  in  a  readable  fashion.  TAPE7 
contains  plotting  data  in  x-y  pairs  (strain  «  x,  stress  *  y  for  strain 
rate  controlled  runs;  time  *=  x,  strain  =  y  for  stress  rate  controlled 
runs) . 

The  format  for  the  input  file  TAPE5  is  as  follows: 

line  1:  Bodner  model,  material  cqeff icients,  free  format,  using  real 
values 

E,  Zq,  Zp  Z2,  m,  A,  r,  n 

line  2:  test  control  parameter,  integer  value  (0,  1  or  2) 

0  implies  strain  rate  controlled  run 

1  implies  stress  rate  controlled  run 

2  implies  run  over 


line  3:  number  of  times  the  controlling  load  (strain  rate  or  stress 
rate)  changes;  integer  value 

line  4:  test  conditions,  free  format,  using  real  values 
Rate,  Time,  Time  step.  Print  interval 

where  rate  =  stress  of  strain  rate  in  this  part  of  the  test 

time  =  time  when  this  test  condition  terminates,  measured 
from  the  start  of  all  loading 
time  step  =  time  increment  used  for  calculations  (for  the 
stress-strain  curves  used  in  this  work,  the 
time  step  was  0.01  sec) 

print  interval  =  time  interval  between  printed  output 
values 

Repeat  line  4  until  there  are  the  same  number  of  test  conditions  speci¬ 
fied  as  there  are  load  changes  indicated  in  line  3  (i.e.  if  line  3  reads 
"3"  there  should  be  three  line  4's  following  it. 

last  line:  test  termination  parameter;  integer  value  of  2,  free  format 

Program  PLOT 

This  program  is  used  to  generate  strain  rate  versus  stress  curves 
(e.g.  Figure  6).  It;  sis  written  in  FORTRAN  5.  Values  for  the  variables 
are  set  by  modifying  the  program  statements  directly.  The  dimensionless 
forms  of  the  coefficients  are  used.  Values  for  m  times  Z^,  n,  AVDO 
(A/Dq) ,  r,  and  Z2VZ1  (Z2/Z^)  must  be  established.  The  form  listed  here 
generates  a  family  of  curves,  one  for  each  different  value  of  Z2VZ1. 

Output  is  into  a  file  called  TAPE5;  x-y  plotting  pairs  for  non- 
dimensionalized  strain  rate  versus  stress  curves  are  written  in  free 


format  onto  TAPE5. 


Program  CYBODR 

This  program  is  used  to  generate  stress-strain  curves  for  cases 
involving  cyclic  loading  (e.g.  Figure  32).  It  is  written  in  FORTRAN  4 
and  is  meant  to  be  run  interactively.  The  Bodner  material  coefficients 
can  be  changed  either  before  compilation  (by  modifying  the  code)  or 
during  execution  when  terminal  inputs  are  solicited.  The  variable  q 
can  only  be  changed  by  modifying  it  in  the  text  of  the  code.  This  is  a 
user-friendly  program  requiring  no  special  instructions  prior  to  opera¬ 
tion.  Output  is  to  TAPE6  in  a  f ormatted-to-read  record  of  time  steps 


and  incremental  variable  values,  and  to  TAPE7  as  free-formatted  x-y 


UUOOUOOUUUUU  UUU  OO  <_>  O  OUUOUUOUU  ouuou  oooo  uuuoou 


PROGRAM  SOONER  (INPUT, OUTPUT, TAPES. TAPES, TAPE7) 


....  BODNERS  MODEL  IS  A  UNIFIED  THEORV  UITH  ONE 

state  uariaile  known  as  the  hardness  parameters... 

....  THIS  CODE  USES  A  FIRST  ORDER  EULER  FORUARD 
integration  scheme  in  order  to  compare  it  to 
OTHER  THEORIES  OF  THE  SAME  TYPE . 


COMMON/A/E , Z0 . 21 . 2 1 . RH. A . R . RN. D0 , UP . EXLHT , EXUPLMT 
_ READ  IN  THE  SOONER  CONSTANTS  AS  SHOUN . 


REU1ND  5 
REU1ND  S 
REUIND  7 


EXLHT • 1 .6  t  EXUPLMT*  70.0 


'  READ<5, *  )  E,20,21,ZI,RH, A,R,RN 
UP-0.0 
SIG-0.0 
EPLAS*0.0 
STRAIN-0.0 
D0-1.E06 
T-0.0 

....UP  IS  THE  PLASTIC  UORK  INITIALLY 
....SIG  IS  THE  STRESS 
....EPLA5  IS  THE  PLASTIC 

- STRAIN  IS  TOTAL  STRAIN  •  ELASTIC  ♦  INELASTIC  COMPONENTS 

- D0  IS  A  SOONER  CONSTANT 

....T  IS  THE  CURRENT  TIME 


Z-Zl-(Zl-ZO)tEXP<-RH*UPl 


- THE  ABOUE  2  EQUATION  IS  USED  TO  DETERMINE  THE  FIRST 

UALUE  OF  2  ONLY. 

URITE(6, 101 )  E,20,Z1,ZI,RH,A,R,RN 

101  FORMAT ( 10X.2X, XE-t, El 1 .4, 2X, 1Z0-X, F6 .2,2X,XZ1-X,F?  .2.3X, IZI -X 
1 ,F£.2,2X,XH-X,F5.2,2X,XA-X,El0.3,2X,XR-X,FS.2.2x,SSMALL  N-X, 
1F4.2,//) 


17  CONTINUE 

N-O^ENOT^STRAIN  CONTROL, KNOUN  STRAIN  RATE  OR  STRAIN/TIME  CURUE 

M-l  DENOTES: STRESS  CONTROL, KNOUN  STRESS/TIME (PIECEWISE  LINEAR)  CURUE 

M-2  INDICATES  A  JUMP  OUT  OF  THE  LOOP  TO  STOP. 

IF(M.EO.O)  GO  TO  IS 
IF (M.EQ. 1 )G0  TO  16 
IFCN.E0.2)  GO  TO  999 

15  CONTINUE 

CALL  EPSLON(T, SIC, STRAIN, Z.EPLAS l 
GO  TO  17 

16  CONTINUE 

CALL  STR£SS<T.SIG.STRAIN, Z.EPLAS) 

GO  TO  1? 

999  CONTINUE 
END 
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FV  V  V  V.V  ^  V 


^ir»4ww?  VT? 


& 


SUBROUTINE  EPSLONCT, SIC, STRAIN. Z, EPLAS ) 
CONMON/A/E.Z0,Z1,ZI.RM,A,R,RN,D0,UP.EXLNT,EXUPLHT 
URITE(6. 150 > 

FORMAT 1 10X, 1X.ITHIS  SECTION  HAS  N-0.FOR  STRAIN  CONTROL , UHERE  THE  S 
1TRAIN/TIP1E  CURUE  IS  IMPOSED.*,//) 

READ(5. * )  N 
UR1TE(6,1®2)  N 

FORMAT<10X,1X,*INPOSED  STRAIN  RATE  CHANGES*.  12, 1X,*TIPIES.*./ ) 

- N  IS  THE  NUMBER  OF  TINES  THE  STRAIN  RATE  CHANCES. 

. ...  NOT  COUNTING  THE  ORIGIN  POINT. 

SSC-S0. 

ESC*®. 5 
URITE(6,103) 

.S22S?UiS2',T5/5Ji.!5ELT5  t*,iox.*stress*,6x,*strain  ratei.sx. 
i*strain*,sx,*plastic*,i2x,*2*;//) 

DO  25  1*1, N 


READC5,*  )  RATE, TINE. DELTAT, PRINT 


...RATE  •  STRAIN  RATE  IN  THIS  PIECEUISE  LINEAR  REGION 
...TINE  *  TINE  WHEN  RATE  FINISHES  AND  CHANCES  TO  THE  NEXT  RATE 
...DELTAT  •  TINE  INCREMENT  THAT  THE  RATE  IS  DIUIDED  INTO 
...PRINT  •  AT  UHAT  TINE  1NTERUALS  A  UALUE  UILL  BE  PRINTED 

P-0.8 

CONST  *  DELTAT*2.0/SQRT(3.0 )*D0 
CONTINUE 

T  -  T  ♦  DELTAT 
P-P+DELTAT 

IF(ABS(SIC).LE.EXLNT)  CO  TO  S . 

X  —  URN+1.  )/(2.*RN> )*(Z/ABS(SIQ))**(2.*RN) 

IF  (ABS(X).CE.EXUPLHT)  CO  TO  5  . 

CONTRB  •  SICN(CONST*EXPCX ),SIG) 

CO  TO  6 
CONTINUE 

CONTRB  *  0.0 
CONTINUE 

EPLAS  ■  EPLAS  ♦  CONTRB 
UPDOT*SIGs( CONTRB J/DELTAT 
SIG*51G+E*<DELTAT»RATE )-E*( CONTRB ) 

STRAIN-SIC/E+EPLAS 

ZD0T-RN*(Z1-Z  )*UPD0T-Z1*A*( (Z-ZI  >/Zl  )**R 
Z*Z+ZDOT*DELTAT 
IF  (P.LT. PRINT)  CO  TO  113 

URITEC6, 104 )  T, DELTAT, SIC, RATE, STRAIN, EPLAS, Z 
URITE(7,» )  STRAIN, SIC 

FORMAT ( 1X,F8.3,5X,F8.4,3X,5E14.4 )  - 

SX1-SIC/SSC 
EPX1-100.ISTRAIN/ESC 
P-0.0 

IF  (T.LT.TINE)  GOTO  10 
T-TIME  -  ‘ 

CONTINUE 

URITE(?,S)  999.999 

RETURN 

END 
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SUBROUT I WE  STRESS ( T. S IG. STRAIN. Z.EPLAS ) 
CONNON/A/E.ZB.Zl.ZI.RN.A.R.RN, D0.UP,EXLHT,EXUPLNT 
URITE16, 151 ) 

151  FORMAT! 10X, IX, (THIS  SECTION  HAS  N-1.F0R  STRESS  CONTROL. UHERE  STRES 
1S/TIME  CURUE  IS  INPOSED.*,// ) 


READC5,*)  N 
URITECb, 112 )  N 

112  FORflAT! 10X, IX, (IMPOSED  STRESS  RATE  CHANGES*. 12. 1X.*TIHES.*,// 1 
SSC*50. 

ESC*®. 5 
URITE(G, 153) 

153  FORMAT! 1SX.*T*.9X, (CELT A  T*, 1®X. (STRESS*. CX. (STRESS  RATE*.6X.*STRA 
1 IN  *.SX,*PLASTlC*.12X.*Z*.'/> 


00  28  1*1  ,N 

READCS.*)  RATE , TINE. DELTAT. PRINT 
CONST •0EITAT*2 . 0/SQRT  <  3 . «  >*00 

11®  CONTINUE 
T.T+OELTAT 
P-P+DELTAT 

IF! ABSCSIG  ) . LE.EXLNT )  GO  TO  165 

X— ! <RN*1.  )/12.*RNJ)*!Z/ABSISIC>>**!2.*RN» 
IF! ABS(X )  .CE.EXUPLMT 1  GO  TO  1®S 
CONTRB*SIGN!CONST*EXP!X  ),SIG > 

GO  TO  108 
10S  CONTINUE 

CONTRB*®.® 

10G  CONTINUE 

EPLAS*EPLAS*CONTRB 

STRAIN*STRAIN+DELTAT*RATE/E+C0NTRB 

UPDOT*SIC*CONTRB/OELTAT 

SIG*E*lSTRAIN-EPtAS) 

ZD0T»RM*!Zl-Z)*UPD0T-2l*A*( IZ-ZI )/Zl)**R 

Z*Z+ZDOT*DELTAT 

PP*P4 • A0OO6001 

IF  (PP.LT. PRINT!  GO  TO  152 

URITEC6. 107 J  T.OELTAT.SIG.RATE.STRAIN.EPLAS.Z 

URITEI7. *  )  T. STRAIN 

10?  FORMAT! 10X, 1X,F8.3,5X,F8.4,3X,5E14.4  > 
SX1-SIG/SSC 
EPX1*100.*STRAIN/ESC 

152  ^CONTINUE 

TTT  *T+0 . 000000001 
IF(TTT.LT.TINE)  CO  TO  110 
124  T-TiriE 
28  CONTINUE 
RETURN 
END 


s 
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PROGRAM  PLOT 


COMMON/UALUE/MZ1 .  AUDO,  2202 1 ,  M.  R.X.  L 

REAL  HZ 1 . AUDO , Z2UZ 1 , N. XP, DUMB V , F . XARRAV, VARRAV. R 

INTEGER  L 

DIMENSION  DUMMY ( 8 ) ,  F  <  2 ) , XARRAV ( S4 ) , VARRAV ( 64 ) 


tuimtmimimtmt 

CREATE  SOME  STARTING  PLACE  FOR  X  IN  THE 
ITERATIVE  SCHEME  TO  FOLLOU 

mummmimtmts 

X  •  .5 

nmmnmttinmtK 

INITIALIZE  NOMINAL  UALUES  OF  THE 
•BACKGROUND*  COEFFICIENTS 


iiiiiitmmmuuittt 

t 

NZl  •  747.859822 
N  *  3.0 

AUDO  •  .000000003 
R  •  7.0 
DO  200  K-1,S 

muimtmmtmmt 
*  INITIALIZE  THE  COEFFICIENT  OF  INTEREST 
mmmmmttmmti 

IF  (K  ,E0.  1)  THEN 
Z2U21  •  0.3 

ELSE  IF  (K  ,E0.  2)  THEN 
Z2UZ1  -  0.35 
ELSE  IF  (K  .EG.  3)  THEN 
Z2UZ1  •  0.4 

ELSE  IF  (K  .EQ.  4)  THEN 
Z2UZ1  •  0.4S 
ELSE  IF  (K  .EQ.  5)  THEN 
Z2U21  ■  O.S 
R  •  7.3 
ENDIF 
S 

V  -  1 .22076312SE-16 
VARRAY(l)  •  V 


CALL  CALC  (Y) 


XARRAV ( 1 ) 


DO  100.1-1,53 
L  *  ‘I  ♦  I 
V  •  2.  x  Y 
CALL  CALC(V) 

xarrayil)  >  x 

VARRAV (L  >  •  V 
100  CONTINUE 

S 

DO  150  ID-1,54 

URITE  (5,0  XARR AY (ID), YARRAY (ID) 
150  CONTINUE 
t 

200  CONTINUE 
t 
t 
t 

999  FORMAT  (FIS. 13) 


SUBROUTINE  CALC(Y) 

COf1HON/UALUE/RZl .  AUD0,  Z2UZ1 , N,R,X,L 

REAL  X.V.nZl,AUDO,Zau21.N,DUIinv,F.XP,R 

INTEGER  M 

DII1ENSION  DUnHY(8).F(2> 


dunnyu  )  •  o.  ini  /  (n  +  i. ) 

do  toe  J-i.i 

ounnY(a)  •  log(e./(Sqrto.  >xv)) 

ALPHA  •  (DUMY(l)  *  DUNNY<2)  )tt(I./(2.*N>) 

Dunnvo)  •  nzi  *  (i. -(alpha  i  x>>  t  x  *  y 

DUrU1Y<4)  •  AUD0  *  ((ALPHA  I  X)  -  Z2UZ1 )  St  R 
ru)  •  Dunnvcs)  -  DunriY<4) 
xp  -  1.1  t  x 


do  se  k-i.so 

n  -  l  -  1 

IF  (N  .CE.  39)  THEN 
X  •  1./  ALPHA 
GOTO  £0 
ELSE 

DUWIY(S)  •  P121  *  (1.  -(ALPHA  I  XP))  t  XP  t  V 

DunnY(£)  •  aud0  r  ((alpha  i  xpi  •  zauzn  tt  r 
F (2 )  -  DUPINY(S)  -  DUHNY(G) 

DUNHY<7 )  •  (F<2 )  *  (X  -  XP))/<F(1)  -  FC2I>  ' 
DUPWY(B)  •  XP  -  DUHNY ( 7 ) 

F(l)  •  F<2) 

X  -  XP 

XP  -  DUNNY(8  > 

IF  (ABSIFd  )  > .  LT.  1  .E-15  )  GOTO  60 


ENDIF 

50  .  CONTINUE 
PRINT  t. 
50  CONTINUE 

100  CONTINUE 
RETURN 
END 


'ROUTINE  DID  NOT  CONVERGE' 
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PROGRAM  CYBODR  ! INPUT, OUTPUT, TAPES. TAPE6.TAPE7 > 


THIS  PROGRAM  SOLVES  BODNERS  LATEST  UNIFIED  THEORY 
FOR  ANISOTROPIC  PLASTIC  FLOU.  Tm£  STATE  VARIABLE 
PARAMETER  2  HAS  BEEN  SPL11  INTO  2  TENSION  AND 
2  COMPRESSION.  BOTH  VALUES  OF  2  ARE  CONTINUOUSLV 
COMPUTED  UHETHER  THE  LOADING  IS  IN  POSITIVE  OR 
NECATIVE  DIRECTION. IRRESPECTIVELY. 

CYCLIC  ANALYSIS  IS  CARRIED  OUT  UITH  TIME  STEP 
AUTOMATICALLY  SELECTED  BY  THE  PROCRAM.  IT  IS 
AN  INTERACTIVE  PROCRAM  FOR  DATA  INPUT.  THE 
NUMBER  OF  CYCLES  ANALYSIS  TO  BE  CARRIED  OUT  MUST  BE 
TEN. 


immnmmmuimmmimmtmmimti 

DEVELOPED  BY  A  .  H  .  RAJENDRAN 
SERVICE  LIFE  MANAGEMENT  SECTION. 

UNIVERSITY  OF  DAYTON  RESEARCH  INSTITUTE. 

mimmnttmuiitimmmmummmmi 


REAL  M.MBAR.N 

COMMON  /  MATCONS/  E ,20 ,21 ,23, 23, M.MBAR.N. A, R.D0.Q 
COMMON  /  LIMITS  /  EXLNT, EXUPLNT, TOLR 

DATA  E/162.SE3/.Z0/1622./.Z1/1795./.Z2/718./, 23/0.0/. 

♦M/.4 l667/.nBAR/e.8/,N/3.0/.A/1.5E-83/,R/7.e/,D0/l.E+e6/ 

Q  •  -.50 
REUIND  6 
REWIND  S 
URITE(S,1001 ) 

URITE  1001 

1001  FORMAT ( ///5X,60( 1H<  J/5X, lHt.SBX, lHt/SX, lHt , HX, t CYCLIC 

♦t STRESS-STRAIN  RESPONSE*. 15X, 1H*/SX,1H* ,58X, 1H*/5X. 1H*,4X, 
♦*CALCULATED  FROM  BODNEfi-PARTON  CONSTITUTIVE  MODELS*. 5X, 1H*/ 
♦5X,1H*.7X.*F0R  ISOTROPIC/ AN ISOTROPIC  MATERIAL  BEHAUI0R*.8X,1M*/ 
♦SX, 1H*,58X,1HI/SX,60( 1H*  )  ) 

CALL  PRINT 
URITE  1003 

1003  FORMAT (///15X,*MATERIAL  CONSTANT  CHANGES  NEEDED?*// 16X. 

♦*(Y/N)  ••>  *) 

20  READ  130.  I 
130  FORNAT(Al) 

I F <  I .EG. *N*  )  GO  TO  250 
IF  <  I  .EO.  *Y*  )G0  TO  50 
URITE  1000 

1000  FORMAT I//1SX,S( 1HS ), t  INUALID  RESPONSE  *,S< 1H* )/18X, 

♦IMLtASE  ENTER  *Y*  >  YESl/aiX.l’H*  •  N0*//B8X,*-->  *) 

CO  TO  20 
50  URITE  1002 

1002  FORMAT t//5X,i7(lH*  ),*  ENTER  MATERIAL  CONSTANTS  I.17UH*)// 
♦8X.IN0TE*  IF  THE  CURRENT  VALUE  OF  A  PARTICULAR  CONSTANT  I St. 
♦/leX.tACCEPTABLE,  ENTER  A  ZERO  IN  THE  FIELD  AND  RETURN*) 

URITE  134. s  . 

134  FORMAT (/&V, 60 ( 1H* )//18X, *£  •  I) 

CALL  INPUT(E) 

URITE  135 

135  FORMAT! 17X, *20  •  *) 

CALL  INPUT<20) 

URITE  136 

136  FORMAT! 17X. *21  •  S) 

CALL  INPUT! 21 ) 

URITE  137 

137  FORMAT! 17X. *22  •  *) 

CALL  1NPUT(Z2) 

URITE  138 

138  FORMAT! 17X, *23  •  *) 

CALL  INPUT1Z3) 

URITE  139 

139  FORMAT! 18X.XM  •  *) 

CALL  INPUT (M ) 

URITE  140 

140  FORMAT! ISX.tNBAR  •  S) 


CALL  INPUT  ( MBAR  ) 

URITE  HI 

HI  f0fir,A?(i8X,iP  •  *) 

CALL  INPUT (A ) 

URITE  H8 

H2  F0RMAT(18X,xR  *  *> 

CALL  INPUT ( R  ) 

URITE  H3 

H3  FORMAT ( 18X,  SQ  •  *) 

CALL  INPUT (0 ) 

£00  CONTINUE 
CALL  PRINT 

VERIFY  CONSTANT  DATA  INPUT 
URITE  H4 

H4  FORMAT ( //SX.XARE  THE  PRINTED  UALUES  FOR  THE  MATERIAL  CONSTANTS!. 

+*  ACCEPTABLE?*'/16X.*(V/N)  *) 

220  READ  130,1 

IF ( I .EQ. *N* )C0  TO  50 
IFd.EO.'Y1  )C0  TO  250 
URITE  1080 
GO  TO  220 

250  EXLnT-1.0  t  EXUPLMT •  70.0  S  T0LR-1.E-1S 


T0LR-1.E-1S 


URITE  1005 

100S  F0RMAT<//5X,22(1H*).*  CYCLE  CONTROL  *.23C1HI >// 1SX. 

♦*ENTER*.?X,**1’  •  stRAIN*/27X,l‘2*  •  5TRESS*//24X.*~>  ») 
READ  X.lPROB 
CALL  INCNTL  (IPROI.E) 

STOP 

END 


SUBROUTINE  ZDOT(ZT,2C,UDRCT.UPD0T,D£LTAT> 

COMMON*  I^MAf CONS/  E. 20,21. 22.23. JJ.MBAR.N, A, R.D0.O 

COMMON  /  LIMITS  /  EXLMt .EXUPLMT, TuLR 

IF(UDRCT.CT.TOLR)  2EFF-2T 

IF (UDRCT.LT.TOLR )  2EFF-2C 

HARD1 «MX (Zl-ZEFF )*UPDOT 

PPPPP-AB5( (2EFF-22J/21 > 

22DCT  •  HARD1-AX21*(PPPPP)**R 

2TDCT  •  (Q+ ( 1 . 0-Q  JXUDRCT )  xZZDOT 

2CD0T  -  ( Q- ( 1 . 0-0  JBUDRCT )X22D0T 

ZT-2T+2TD0TXDELTAT 

2C-ZC+2CD0TXDELTAT 

RETURN 

10  URITE(6,1000)CH1,CH2 

U9ITE  1000, CHI, CH2  _  ...  ,  ,.u 

1000  FORNAT<//SX.22( 1M* ),*  END  OF  PROGRAM  I ,22C 1H» )//10X, 

+ JTHE  EXPONENT  OF  Th£  RECOVERY  TERMS  HAVE  BECOMEX/10X. 

♦  XHEGATIUE  IN  SUBROUTINE  ,ZD0T,*//16X,*CM1  *  B.E15.4/16X. 

♦  XCH2  -  B.E1S.4,//)  " 

STOP 

END 


SUBROUTINE  INPUT  <  CNST  )  . 

READ  *,  CNSTMAT 

IFU.hSTnAT.LE,  0.000001)  RETURN 

CNST  ■  CNSTMAT 

RETURN 

END 


SUBROUTINE  PRINT 

REAL  n  .MBAR.N  „ 

COMMON  /  MATCONS  /  E.ZO,Zl.Z2,23,n.MBAR,N,A.R,D0.O 
COMMON  x  LIMITS  /  EXLMT, EXGcNMT , TOLR 
URITE  iee,E,2O,Zl,Z2,Z3,N.niAR.N,A,R,D0.O 
URITE (£,  100  )£,* 0,21 ,22, Z3,M, MBAR, N,A,R, D0, Q 
100  F0RMAT<//SX,i5<iAb>,*  6ur6ENT  MAtE^IAL  CONSTANT  VALUES  *,  .. 

♦13<1H*>///1$X.*E  •  *,F1S.4/17X.*Z0  •  JF15.4/17X.0Z1  ;*.flS.4/ 
♦17X.IZ2  •  *.Fis.4/17fc,*Z3  •  I.F1S.4/18X.BM  •  «,F6.3/i5X.*MBAR  *  *. 
♦F£.J/18X,tN  •  X,Ffi.3/18X,IA  •  f ,Ffi.4/18X,*R  •  *,F£.3/1?X,«D0  •  *, 
♦F7.1/I8X.B0  •  B.F6.3/) 
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RETURN 

END 


SUBROUTINE  INCNTL  <  IPROB.E  ) 

URITE  100 

100  FORMAT (//1SX.XT0TAL  NUMBER  OF  CYCLESS//B4X,X-->  X) 

READ  t , NCVCLE 
WRITE  120 

120  F0RNAT(//1SX.IRAXIHUM  AND  RINIRUR  STRAIN  OR  STRESSX//24X.*— >  X) 
READ  t,  SNAX  ,  SHIN 
URITE  140 

140  FORMAT (//15X.SNURBER  OF  CVCLES  PER  RINUTEX//24X,X-->  t> 

READ  *,CPR 
CP  S-CPR/60.0 
PERIOD  *  1.0  /  CPS 

RATE  •  <  SC1AX  ♦  SHIN  )  /  PERIOD  X  2.0 

TOTCV  •  NCVCLE 

TMAX  •  PERIOD  *  TOTCV 

IF ( IPROB.EO.2)  TSTEP-0.000HE/RATE 

IF(lPROB.EO.l)  TSTEP-0. 0001/RATE 


CALL  STRSTN  < I PROS, ZT.ZC.R ATE, TSTEP.TNAX, PERIOD. SNAX. SNIN ) 


RETURN 

END 


SUBROUTINE  RATECHG  <  RATE, T.TSTEP, PERIOD, TCHECK.IDENT  , DTOLD 
♦  .SriAX. SHIN. PCHECK) 

TOLR-1.0E-1S 

IF (T . LE.TOLR )  RETURN 

TTT-TCHECK-TOLR 

TERP-ABS(RATE) 

IF ( IDENT.EO. 1 )  PCHECK -SRAX/TERP 
IF ( IDENT .EQ.2)  PCHECK-SRAX/TEHP 
IF< IDENT.EO. 3)  PCHECK-SRIN/TERP 
IF( IDENT.EQ.4 )  PCHECK-SNIN/TERP 
I F ( TTT . LT . PCHECK )  RETURN 
TSTEP-DTOLD 
TCHECK  •  0.0 
URITEIG.llli) 

1111  FORRATtSx, t  999  999  999  999  999  999  999  X.) 

OO  TO  I  1,2. 3. 4  )  ,  IDENT 

1  RATE— RATE 
IDENT-2 
RETURN 

2  I  DENT-3 
RETURN 

3  RATE— RATE 
I DENT -4 
RETURN 

4  I DENT -1 
RETURN 
END 


SUBROUTINE  STRSTN  (  IPROB,ZT,ZC,RATE,TSTEP,THAX.PERIOD,SIW<,SniN  ) 
REAL  R.RBAR.N 

CORRON  /  RATCONS/  E,Z0.Zl,Z2,Z3, R.RBAR.N, A, R.D0.Q 
CORRON  /  LIMITS  /  EXLRT , EXUPuRT , TOLR 

INITIALIZE  THE  VARIABLES 

ZT-Z0  >  ZC-ZB  B  DTOLD-TSTEP  •  ISTOP-1  t  EPINCR-0.0 
UP-  B.B  B  SIG  -  B.B  S  EPLAST  •  B.B  B  ETOT  •  0.0 
T  •  0.0  B  IDENT  •  1  B  TCHECK  •  0.0  S  TPRIMT  •  0.0 
PCHECK • SR AX/RATE 

PRINT  HEADING 

URITE  (6.200) 

200  F ORRAT ( ///5X , 9S ( 1 Ht ) // 

♦IBX.XTIREX.SX.XTIRE  STEPB.SX.XSTRAIN*. 1BX.XSTRESS*. 

♦SX.XINELAStIC  STRAINS. 3X,*Z-tENSI0NI,5X.XZ-C0nPRESSI0NX/' 
♦SX.95(1HS>//) 
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PRINT •  PERIOD  /200.0 
10  CONTINUE 

CALL  RATECHC  (  RATE, T, TSTEP, PERI OD.TCHECK, I DENT, DTOLD, 
♦SNAX, SHIN. PCHECK > 

CALL  TINECHC l EPINCR, IDENT. TSTEP, ISTOP  ) 

CALL  T1N£ADJ< PERIOD, TSTEP, TCHECK, DTOLD, IADJ.PCHECK ) 

CONST*  TSTEP  12./  SORT  (  3.0  )*  DO 

T-T+TSTEP 

TCHECK-TCHECK+TSTEP 

TPRINT-TPRINT+TSTEP 

IF (AISCSIC ) . LE.EXLHT )  CO  TO  100 

IF (UDRCT.CT.TOLR )  Z-ZT 

IF ( UDRCT . LT . TOLR )  Z*ZC 

X—  ((NM.  )/(2.*N>)*(  ((Z/ABS(SIG))**2)**N> 

IF ( ABS(X  )  .CE .EXUPLNT )  CO  TO  100 
EPINCR 'SIGN (CONST  *  EXP(X).SIC) 

CO  TO  110 
100  CONTINUE 
EPINCR-0.0 
110  CONTINUE 

EPLAST-EPLAST+EPINCR 
WPDOT*SIC*EPINCR/TSTEP 
CO  TO  (120,130).  1PR0I 
120  SIC*SIG+EI(TST£P*RATE-EPINCR) 

ET0T-SIC/E4EPLAST 
CO  TO  140 

130  ETOT-ETOT*TSTEP*RATE/E+EPINCR 
SIG-EKETOT-EPLAST) 

HO  CONTINUE 

UDRCT*SIG/ABS(SIC) 

CALL  ZDOT( 2T . ZC, UDRCT, UPDOT, TSTEP ) 

TTT-TPRINT+TOLR 
IFtlMDJ.EO.l  )  CO  TO  H5 
1F(TTT.LT. PRINT)  CO  TO  ISO 
MS  CONTINUE 

WRITE (7, 222  )  T. TSTEP. ETOT, SIC, EPLAST.ZT.ZC 
WRITE  (S,*)  ETOT,  SIC 
222  F0RnAT(5X,7(E12.5,4X ) ) 

TPRINT  «  0.0 
ISO  CONTINUE 
TTT-T+TOLR 

IF(TTT.LT.TNAX)  CO  TO  10 

RETURN 

END 

SUBROUTINE  TINECHG( A, IDENT, TSTEP, ISTOP ) 

TOLR-l.OE-06 
CR*ABS( A ) 

CO  TO  (10.20.10,20)  ,  IDENT 
10  IF (CR  .LE.  TOLR)  RETURN 
IF( ISTOP  .EQ.  2)  RETURN 
TSTEP* TSTEP/ 10.0 
ISTOP-2 
RETURN 

20  IF(CR.CE.TOLR)  RETURN  " 

IF (ISTOP  .EO.  3)  RETURN 

TSTEP ‘TSTEP* 10.0 

ISTOP-3 

RETURN 

END  -  s 


SUBROUTINE  TIHEADJ(PERIOD, TSTEP, TCHECIC.DTOLD, IADJ.PCHECK  > 

DT 1  •  PCHECIC-TCHECK 

IADJ*0 

IF (DTI  .CE.  TSTEP)  RETURN 

IADJ-1 

DTOLD-TSTEP 

TSTEP-DTl 

RETURN 

END 
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